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The study of production of quark-gluon plasma (QGP) in the laboratory is fascinating for a 
variety of reasons. It is interesting beacuse, this phase did exist for first few microseconds in 
the early universe and perhaps also exists in the stellar matter at high density. Further, a study 
of QGP provides a test of QCD in a new scale (high energy density) in its non-perturbative 
regime. It is also expected to improve our understanding of other fundamental concepts such as 
the nature of confinement and and chiral symmetry. 

From quantum chromodynamics as input dynamics for strong interaction between quarks and 
gluons, we can directly predict this state of matter by large scale computer simulations. The 
results of computer simulations provided a quantitative confirmation of the idea that hadronic 
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matter at sufficiently high density will undergo a transition to a state of matter of deconfined 
quarks and gluons, i.e. quark-gluon plasma On the experimental side, there has been much 
effort in the direction of detecting this phase of matter. In the last decade two experiments at 
SPS of CERN and AGS of BNL were done which succeded in achieving a high energy density, 
although the creation of QGP itself is not confirmed. Future experiments at relativistic heavy- 
lon collidres (RHIC) of BNL and large hadron colliders (LHC) of CERN will hopefully able to 
produce energy density high enough for the production (and equilibration) of QGP 

The most difficult task in these experiments is the detection of quark-gluon plasma. This is 
because, the plasma exists for a very short time in a very small volume and a direct detection is 
not possible The signatures suggested are necessarily indirect and the prominent ones are. (1) 
J/0 suppression, (2) electromagnetic probes such as dilepton and direct photon production, and 
(3) strangeness enhancement. The measurement of these signatures involve both theoretical and 
experimental constraints. The main experimental difficulty lies in singling out the contributions 
arising from the non-qgp sources. The major theoretical difficulty in the calculation of these 
signatures arise from the poorly understood space time evolution of the quark-gluon plasma in 
these ultra relativistic heavy-ion collisions. 

The space time evolution of the quark-gluon plasma evolution at these experiments hats been 
generally described by hydrodynamics assuming that the equilibration is almost instantaneously 
reached. In reality, the various stages by which the complete evolution of quark-gluon plasma 
are to be described are, i) pre-equilibrium, ii) equilibrium, where one actually studies the 
thermalised quark-gluon plasma, iii) cooling and iv) hadronisation. It is crucial and very much 
necessary to study the pre-equilibrium stage in the ultra relativistic heavy-ion collisions, for it has 
direct impact on the plasma evolution and also on the measured signatures in the experiments. 
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In this thesis we study the pre-equilibrium evolution of quark-gluon plasma in URHIC by 
employing color flux-tube model. 

In chapter 1 we introduce the topic and provide the necessary background It contains an 
overview of the theoretical and experimental status of the field. A brief description of the 
proposed signatures is presented in this chapter. 

chapter 2 contains an extensive description of the pre-equilibrium evolution of quark-gluon 
plasma within color fiu.x-tube model. In this model the two nuclei that undergo a central 
collision at RHIC and LHC are highly Lorentz contracted as thin plates. These two Lorentz 
contracted nuclei pass through each other and, in the process, acquire a non-zero color charge 
{{Q) = 0, (Q^) 0), by exchanging soft gluons. This leads to creation of a chromoelectric field 

between two receeding nuclei, which polarise the QCD vacuum and produces quarks and gluons 
via Schwinger mechanism (with suitable adaptions and modifications) The quarks and gluons 
so produced collide with each other and get accelerated by the background chromoelectric field. 
In this chapter we set up the transport equations in an extended phase space (which include 
color in addition to coordinate and momenta) to study the production and equilibration of 
quark-gluon pla.sma m ultra relativistic heavy-ion collisions. 

In chapter 3 we solve the non-Abelian transport equations for qq plasma in SUf2) gauge group. 
In this exercise which is essentially a warm up for the more realistic study which we perform in 
the next chapter, it is shown that the non-Abelian features are crucial and the so called Abelian 
dominance which was studied so far is not justified. We obtain the distribution function of 
quarks and gluons which carry explicit color charges those are abscent in any other study. The 
physical quantities such as local energy density, number density and temperature of the plasma 
are determined for this qq plasma. 
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In chapter 4 we again solve non-Abelian transport equations to include gluons in addition to 
qq in SU(3) group which is the realistic gauge group for strong interaction We made a detailed 
comparison with the SU(2) results as well as with other p-QCD based approaches such as parton 
cascade model. 

A calculation of the dilepton rate is presented in chapter 5 using the distribution of quarks 
and antiquarks within color flux-tube model. It is shown that the pre-equilibrium production of 
dilepton pair is dominant for very small transverse mementum of the dilepton pair, whereas the 
Drell-Yan contribution is dominant over large transverse momentum of the dilepton pair 

In chapter 6, we propose a new source term for quarks and gluons as is appropriate for URHIC. 
We treat the time dependence of the background field dynamically, and also project the amp- 
litudes to a coherent basis in the extended phase space, so as to obtain a legitmate classical 
source term. We show that the ‘source’ is also necessarily a sink both of which exhibit special 
characters in the phase space. We make a detailed comparison with the source term which are 
obtained using rather ad-hoc prescriptions. 

In chapter 7 we have studied J/tv suppressions in a thermally equilibrating quark-gluon plasma 
at RHIC using short-distance QCD. The evolution of the system is taken from parton cascade 
model. The suppression is found to be 100 percent at RHIC energy. A discussion on the validity 
of the models, which describe the evolution of the quark-gluon plasma, is presented from this 
study. 

chapter 8 contains a very brief literature survey of the other models available for the study of 
the pre-equilibrium evolution of quark-gluon plasma. It also contains a very brief description 
of the lattice approaches towards this field. 
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Chapter 1 


Introduction 


In recent times particle physicists and nuclear physicists have evinced a lot of interest in studying 
a new phase of matter, namely, quark-gluon plasma(QGP). It is theoretically established that 
this state of matter consists of unbound quarks and gluons at a temperature about 200 MeV. 
Such a phase did exist in the early stage of the evolution of our universe, ~ 10““^ second after 
Big Bang. As the universe expanded and cooled down, it went through a phase transition to 
the hadronic matter that we see today. To study this phase and to re-create the early universe 
in the laboratory is a challenging task for both the theorists and the experimentalists. 

Quantum chromodynamics (QCD), which successfully describes the strong interaction between 
quarks and gluons, has predicted the possibility of such a phase. Within QCD, quarks carry 
three types of color charges which interact among themselves via gluons, which are the non- 
abelian gauge bosons. Unlike the Abelian theory where photons do not carry any electric charge, 
gluons in the non- Abelian theory (described by the gauge group SU(3)) carry a charge which 
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IS a vector m eight-dimensional space. Within this theory, quarks belong to the fundamental 
representation of SU(3) group and gluons to its adjoint representation This theory explains 
many phenomena observed at the present day experiments. 

However the full range of QCD is not solved till today. This is because, the perturbative 
calculation is not applicable over the large distance scale due to the nature of the strong in- 
teraction between quarks and gluons. This regime needs to be studied within the frame work 
of non-perturbative QCD So far, computer simulations on the lattice[l] provide the only reli- 
able method available for the non-perturbative study. More importantly lattice QCD predicts a 
phase transition between the hadronic state to a quark-gluon phase. It also predicts that at a 
critical temperature Tc the dimensionless energy density measure, e/T'^, increases abruptly from 
a value near that of an ideal pion gas (~ 1) to a value that of an ideal quark-gluon gas. This 
increase in the energy density, Ae, represents some thing like “latent heat of deconfinement”. 
This suggests that for a stable and large system of vanishing net baryon density, there occurs 
a deconfined phase of the quarks and gluons, called quark-gluon plasma. This deconfinement 
phase transition is expected to occur at a temperature of the order Tc = l50MeV corresponding 
to an energy density of Cc = 1-3 GeV//m^. Energy densities of this magnitude which did exist 
in the early universe is expected to be obtained in small volumes (~ 100/m^) and for a short 
time (~ 5 — 10/m) in high energy heavy-ion collisions. Abundant multiple scattering among the 
incident nucleons in the heavy-ion collisions and the many secondaries produced during these 
processes are expected to give rise such an energy density. 

As no direct observation exists in these experiments for the initial energy density released, it 
has to be determined from the distribution of the final state hadrons. The formula obtained 
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from the Bjorken model[2] which is generally used to estimated the energy density is 

_ dErldy 

7tB?To 


( 1 . 1 ) 


where tq = 1/m being the formation time and R = I.ISA^^^. The quantity dEr/dy in the 
above equation is proportional to the final multiple distribution dN/dy which increases as the 
number of binary nucleon-nucleon interaction increases. In order to obtain such an energy 
density through large multiple collisions, it is required that the two nuclei in the collision should 
be large and accelerated to high energy. 


During the last decade such experiments involving heavy-ion collisions have been done at the 
AGS[3] of BNL and the SPS[4] of CERN with Eiab = 14 and 200 A GeV. Both are fixed target 
experiments and hence produce a low centre of mass energy {Ecm)- Although the results of 
these experiments suggest an extremely rich physics which can not be explained by a straight 
forward extrapolation of nucleon-nucleon collisions, it is still not clear whether (an equilibrated) 
quark-gluon plasma is created in these experiments, in particular at SPS (Pb-Pb collisions at at 
Ecm = 17 GeV). The future heavy-ion collider experiments at RHIG (.Au-Au collision) of BNL 
and LHC (Pb-Pb collision) of CERN with Ecm = 200 A GeVand 5.5 A TeV respectively will 
provide the best opportunity to study such a phase. These values of high energy imply that 
at the moment of maximum nuclear overlap the energy densities are of the order of e, ~ 250 
(7500) GeVj frm? at RHIC (LHC), which is much higher than the critical energy density for 
deconfinement caluculated from lattice QCD. 


The most important task on hand is a detection of QGP, assuming that it is produced. As 
the plasma exists for a short time (~ 5 — 10/m) over a small volume (~ 100/m^), it is not 
directly accesible in these experiments. Indirect signatures are proposed to detect this phase 
and the promments among them are ; (1) J/V’ suppression[5], (2) electromagnetic probes such 
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as dilepton and direct photon production[6, 7], and (3) strangeness enhancement [8]. 

• J/4> suppression 

According to the proposal of Matsui and Satz, the J jtj} production is suppressed m the presence 
of a deconfined QGP medium due to Debye screening. At high temperature the Debye screening 
length which is calculated from lattice QCD is smaller than the J/i/j bound state radius. This 
forbids the binding of cc to J/'^. Although suppression of J /ijj is seen in heavy-ion collisions, 
still it is not clear that it is due to screening only This is because, Jjijj is also supressed 
in p-A collisions and A-B collisions at SPS, where there is no QGP phase. This can be seen 
from Fig- 1.1 where the survival probability of J/i/) is plotted against various nuclei in different 
experiments. 

The suppression of Jjxp m the above cases is due to the presence of a nuclear medium. The 
inelastic collision between cc and nucleon explains the data in p-A collisions and A-B collisions 
at SPS. Although a recent measurement of J/V’ by NASO experiment (see Fig-1.1) is not ex- 
plained by the presence of a nuclear medium, it is still not clear whether this supression can be 
understood in terms of the presence of a deconfined medium. Recently Kharzeev and Satz have 
employed the results of Bhanot and Peskin (1979) to study the interaction of heavy quarkonium 
with gluons. This work provides a dynamical ba^is to study the J/i/> suppression by using 
short-distance QGD. Within this formalism J/'^ dissociates to DD due to the interaction with 
hard gluons which are present in the deconfined medium. This mechanism is expected to give 
rise a huge suppression of J/i/) at RHIC and LHC[9, 10]. Besides this, various authors have 
devloped different proposals for explaining the NASO data without assuming any deconfined 
medium[ll-14]. 
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1/3 1/3 

A +B 

Figure 1.1; The Jjif) survival probability after absorption through nuclear matter. Here A and 
B are the mass number of the colliding objects. The full line is the sutvival probability for 
the proton-nucleus (nucleus-nucleus) systems, calculated with a cross section a = 6.2 mb. The 
marks are the experimental datas, as mentioned in the figure. 


In short, although a supression of J/i/? is seen in these existing experiments, no clear conclusive 
evidence for the existence of quark-gluon plasma can be assured. A detailed description of the 
recent theoritical and experimental developements on production and its suppression is 
presented in chapter 7. 


• Electromagnetic probes 

Dileptons and single photons have long been proposed as useful probes of the plcLsma[7]. Once 
produced, they hardly interact with the strong matter and thus carry the details of the circum- 
stances of their production. The major processes for dilepton and photon production in URHIC 
are, (i) hard parton scatiering, (ii) electromegnetic decay of hadrons and (iii) production from 
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partons present in QGP (thermal production). The hard parton scatterings produce high pt 
lepton pairs and photons which are calculated from pQCD. The leading order processes are 
quark-antiquark anihillation and quark-gluon compton scattering. This is simillar to hadronic 
collisions. On the other hand, electromagnetic decay of produced hadrons is the main source 
of dilepton and photon production which obscure the signal of interest (thermal production), 
which are produced from QGP. More importantly, the thermal emission of dilepton and direct 
photon IS not well understood either theoretically (due to lack of proper description of the space 
time evolution of QGP) or experimentally. These thermal lepton pairs and thermal photons 
are expected to be important signatures of QGP because, their measurement is probably the 
only possibility to have a direct access to the transient dense and high temperature phase of a 
heavy-ion collision. Using thermal parton distributions, the emission rate of thermal dilepton 
and photon is calculated in a simillar way as the high pt direct photon and dilepton calcula- 
tion using the structure functions of nucleons inside nucleus. It is expected that in high energy 
heavy-ion collisions such as at RHIC and LHC, the thermal production will be more than the 
other processes. 

At present, a lot of effort are being made both theoretically and experimentally to understand 
the thermal production of dileptons and photons from the quark-gluon plasma. This is one of 
the hot topic in the search of QGP. A detailed description of this topic is presented in chapter 
5. 

• Strangeness enhancement 

Strangeness enhancement was proposed as a possible signature of the QGP state, almost 17 
years ago[15]. Recall that QGP is formed when the matter density is 10 times greater than 
that of the nuclear matter. The energy corresopnding to the temperatute which is greater than 
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that of the critical temperature is larger than the mass of the ss pair (~ 300 MeV). Therefore, 
the ss will occur at almost the same rate as uu or dd in the QGP at high temperature. In the 
hadron gas the strange particles are produced through associated production during the collision 
of two typically non-strange hadrons The threshold energy to produce this strange particle is 
determined by the mass of the strange hadron pair. The reaction with the lowest threshold is 
p + n— > A° + A'’"'' + n which requires 671 MeV. On the otherhand. as we have discussed, in 
a QGP the threshold for the strangeness production is only around 300 MeV. Therefore, the 
strangeness production in the QGP is favoured from the production in the hadron ga5. The 
strangeness enhancement has been investigated by many experimental groups[16] and recently 
it has been one of the important topics in relativistic heavy-ion collisions. 

The signals discussed above had been initially studied by assuming an equilibrated QGP. 
However, the measurement of these signals in the actual experiment, involve the contributions 
from all the stages of evolution of the system. In order to make a convincing case for observation 
of QGP, the observation of several signals must be supported with a reasonable description of 
the interactions of the final state particles. The formation and study of quark-gluon plasma 
in future experiments, will crucially depend upon a detailed microscopic understanding of the 
plasma. Questions relating to its formation, evolution and hadronisation deserve a much better 
study. 

Of particular importance to us is the production and the evolution of quark-guon plasma after 
two nuclei collide with each other. It is e.xpected that in the center-of-mass frame, two highly 
Lorentz contracted nuclei pass through each other and deposit their kinetic energy gradually 
behind them[17]. A qualitatively accepted space-time evolution of these reactions is schematic- 
ally illustrated in Fig-1.2. The various stages by which the complete evolution of quark-gluon 
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Figure 1.2: A typically accepted space time evolution picture of the nucleus-nucleus collisions 
in ultra relativistic heavy-ion collisions 


plasma is described in URHIC are, i) pre-equilibrium, ii) equilibrium, where one actually stud- 
ies the thermalised quark-gluon plasma, iii) cooling and iv) hadronisation. In the very early 
stage of this collision, the partons are being produced and they collide with each other in the 
pre-equilibrium stage to form the equilibrated quark-gluon plasma in the stage (ii). It is expec- 
ted that most of the entropy and transverse energy are produced in this satge. The dynamics 
of partons during the first few fermi determines the initial (and boundary) conditions for the 
hydrodynamic expansion of the plasma phase. The validity of the calculations invoh'ing hy- 
drodynamic expansion, assuming very fast local thermodynamic equilibrium after the partons 
are produced, is to be checked and a proper study of pre-equilibrium dynamics is needed. In 
particular, we want to know how the partons evolve, under what conditions they thermalise. 
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the time scale of equilibration and the typical energy density soon after the thermalisation 
Without these inputs there are always uncertainties in choosing the equilibration time and the 
initial temperature which crucially determines the predictions of all the signatures. To study 
these questions systematically in a realistic model, a detailed study of the pre-equilibrium stage 
during the very early stage of the high energy nuclear collisions is necessary. 

As mentioned, the production of an equilibrated quark-gluon plasma in stage (li) crucially 
depends on the pre-equilibrium evolution, and hence all the signatures m the experiments. 
Indeed it has been pointed out[9. 10] that J/?/i suppression which is generally attributed to the 
equilibrium stage can also arise from the pre-equilibrium stage, an observation which is true for 
strangeness enhancement as well. The contribution of the pre-equilibrium stage to the dilepton 
production and the direct photon production would be even more pronounced, and a careful 
study of the above signals ought to shed light on this stage of the plasma. 

In this thesis we discuss the production and the equilibration of quark-gluon plasma(QGP) in 
URHIC by employing the Flux-tube Model[18, 19] which is a generalization of the familiar Lund 
string model[20] widely used for e'^e“ and p — p collisions. A chapter wise summary is given in 
the synopsis (see page v). 



Chapter 2 


Pre-equilibrium evolution of quark-gluon 
plasma 


2.1 Introduction 


The space-time evolution of quark-gluon plasma created in ultra relativistic heavy-ion collisions 
is governed by hydrodynamics, once it reaches equilibrium The relativistic hydrodynamic 
equations are; 


and 


= 0 

(2.1) 

i9^n^(r) = 0, 

(2.2) 
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where T^‘'{r) is the local energy-momentum tensor and n^(r) is related to the local number 
density n(r) by n(r) = n^{r)u^{r) with u^(r) being the flow velocity. In the absence of any 
dissipation the energy-momentum tensor is written in terms of energy density e(r) and the 
pressure p{r) as 

T^‘'(r) = -p^'"p(r) + (e(r) + p(r))u^(r)u‘"(r). (2.3) 


where e(r) = T^^u^[r)u''[r). If one further neglects the transverse expansion in the early stage 
and considers only the longitudinal expansion, we then arrive at the equation[2]: 


dc e -h p 
dr T 


= 0 , 


(2.4) 


first proposed by Bjorken to describe space-time evolution of QGP In the above equation, 
r = ^(t^ — z'^) is the proper time. 


As these hydrodynamical equations are applicable only in equilibrium, which is attained at 
certain time after the two nuclei have collided with each other it is desirable to determine 
the equilibration time by some microscopic approach. After all, the initial conditions that 
govern the hydrodynamical evolution (Eqn 2.4) can only be obtained from a careful study of 
the pre-equilibrium regeime in heavy-ion collisions For an ideal quark-gluon plasma of ma.ssless 
partons, the equation of state is given by[17] 


p = a{T)T^ 


and 

e ~ 3a(r)r^ (2.6) 

s = ~ Aa{T)T^ (2.7) 

and 


n = b{T)T^ 


(2.S) 
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Here s is the entropy density, T is the temperature, n is the number density of the plasma and 
the quantities are given by 


TT^ 21 

«(^) = ^(16 + yn/) 

(2.9) 

Kr)= ^(16 + 9n,) 

TT 

(2.10) 


with rif being the number of flavours. It follows from the above equation of state together with 
Bjorken’s scaling picture that the temperature of the plasma scales like T{t) cc energy 

density like e(r) a and number density like n(T) cc r~h 

The major difficult task is the detailed study of the poorly understood pre-equilibrium stage of 
the collisions Crudely speaking, this pre-equilibrium stage is the duration between the moment 
when the two nuclei just pass through each other to the actual establishment of an equilibrated 
quark-gluon plasma. While most of the earlier work refers to the equilibrium stage studied by 
employing the hydrodynamic evolution, it is necessary to study the pre-equilibrium stage for 
the above mentioned reasons. Indeed, without a proper study of the pre-equilibrium stage there 
are always uncertainties in the predictions of all the signatures in the ultra relativistic heavy- 
lon collisions (at RHIC and LHC). Any pre-equilibrium study also presupposes a production 
mechanism for QGP, and here, we study the pre-equilibrium evolution of quark-gluon plasma 
in ultra relativistic heav\-ion collisions (URHIC), within color flux-tube model. The color flux- 
tube model is a generalization of the familiar Lund string model widely used for e'^e' and p — p 
colhsions[20]. .A brief description of the model is given below. 
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2.2 Color Flux- tube model 

In this model the two nuclei that undergo a central collision at RHIC and LHC are highly lorentz 
contracted as thin plates[21, 22] (The two nuclei m collision at RHIC and LHC travel almost 
at a speed that of light). These two lorentz contracted nuclei pass through each other and, in 
the process, acquire a non-zero color charge {{Q) = 0, (Q^) ^ 0), by exchanging soft gluons. 
So one may figuratively call such nuclei after collision as color capacitor plates between which a 
strong chromoelectnc field is created. The nuclei which act eis color capacitor plates produce a 
chromo-electric field between them [21, 22] The strength of the field which naturally depends on 
the strength of the color charge residing on the plates cannot be fixed from first principles. We 
can only fix that phenomenologically, say by identifying the field energy with the energy in the 
central region as estimated by Bjorken[2] This strong electric field creates qq and gluon pairs 
via the Schwinger mechanism which enforces the instability of the vacuum in the presence of 
an external field. The partons so produced, collide with each other and also get accelerated by 
the background field. The mutual collisions drive the system towards equilibrium with suitable 
modulation from the background acceleration Within this model a set of transport equations 
are solved to study production and equilibration of quark-gluon plasma in URHIC. A brief 
description of works done on this model is given below. 

In the approach of Baym[lS], the partons are assumed to be created soon after the two nuclei 
have collided with each other, and are further assumed to equlibrate almost instantaneously 
thereafter (with a very small collision time). This quark-gluon gas then proceeds to expand hy- 
drodynamically. Subsequent work of Kajantie and Matsui[19], incorporated a dynamic particle 
production by introducing a source term via the Schwinger mechanism[23] in the Boltzmann 
equation, and they also studied the equilbration in the relaxation time approximation. On the 
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other hand, Bialas et a/[24] generalized the Baym analysis by introducing the effect of the back- 
ground field on the otherwise hydrodynamic flow. A self consistent study of the system was 
carried out by Banerjee et a/[25] who combined the eff'ect of both the background field and the 
collisions on the quark-antiquark system. Finally, by employing the same analysis as of Banerjee 
et al, Asakawa and Matsui[2S] have studied the rate of dilepton production by qq annnihilation. 
as a function of (proper) time .A.11 the above analyses ignore the gluon component. 

However some of the crucial features which have not been considered in these earlier studies 
remains to be justified. First of all in the analyses of Kajantie et a/[19] and Banerjee et a/[25], the 
source term acquires a time dependence by virtue of the time dependence that the electric field 
suffers because of the particle production The production of particles is then strictly go\erned 
by a time dependent field, in which case a perturbative mechanism will take over from the non- 
perturbative Schwinger mechanism. Note that the latter holds only for a uniform constant field. 
All the previous studies employing the flux tube model employ the non-perturbative expression. 
2) The source term is derived for QED processes and needs to be redone for QCD processes at 
hand It is apparent that the QCD effect will show up most manifestly in the gluon production; 
there is no corresponding counterpart in QED. In other words it is necessary to obtain a source 
term for gluons. 3) It is well known from earlier studies of Yang-Mills equations that there 
are inequivalent gauge field configurations w'hich yield the same field tensor[26], and hence the 
same energy momentum tensor. This feature which is again absent in the Maxwell case has to 
be taken into account in this study. Recall that the lattice studies at any given energy density 
implicitly sum over the contributions coming from all such inequivalent configurations[27]. 4) 
Finally, the inherent non-linearity and the local gauge invariance which characterize the Yang- 
Mills equations lead to an added complexity; the equations do not possess a unique solution, 
even after gauge fixing, and even if one considers only static configurations. Consequently, one 
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needs to scan, with some suitable weightage, all the solutions for the given initial conditions, 
and sum over all such configurations It is not clear if any procedure is known that allows an 
implementation of this requirement 

We have attempted here to fill this gap partially within color flux-tube model by solving transport 
equations m an extended phase space 


2.3 Transport equation in an extended phase space 


In the study of production and equilibration of quark-gluon plasma in ultra relativistic heavy-ion 
collisions one needs to know the distribution of quarks and gluons, f{x,p), in the usual phase 
space. The determination of all physical quantities in the experiments crucially depend on this 
distribution functions of quarks and gluons. At equilibrium the above distribution functions 
are the standard Fermi-Dirac distribution function for the quarks and anti-quarks and Bose- 
Einstein distribution function for gluons. In the non-equilibrium stage the time devlopement of 
the phase-space distributions fa{x,p) for the partons of species a = q/^qj, g (f being the quark 
flavors) is described by the relativistic transport equation (Boltzmann equation); 


-k Ff,^p''d^\ /a(x,p) = Caix,p) + Sa{x,p) 


(2.1i) 


with some known initial condition of fa{x,p) at t = to, when the two nuclei begin to overlap at 
RHIC and- LHC. The first term in the left hand side of the above equation corresponds to the 
usual convective flow and the second term is the Lorentz force term. 5 and C on the right hand 
side correspond to the source term for parton production and collision term among the partons 
(described in chapter -3). 
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The above equation is written in the usual six dimensional phase space and does not carry 
any information on the time evolution of the color charge, which is a dynamical quantity for 
non-Abelian plasma. Thus, as such, these transport equations which are valid for an abelian 
plasma are not applicable to quark-gluon plasma. It is our purpose to emphasize and study a 
crucial feature of the quark-gluon phase, tinr, the inherently non- Abelian nature of the plasma. 
(Note that all the earlier macroscopic approaches [18. 19, 24, 25, 28] are in the so called Abelian 
approximation). Put simply, the non-Abelian features are completely ignored in the space i.ime 
evolution, and the Boltzmann equation describes the dynamics which is essentially that of a 
Maxwell fluid, be it in the source term or in the background field term. The approach of Geiger 
and Kapusta[17] does incorporate the non-Abelian features in their basic Feynman diagrams, 
but IS neverthless incomplete in that the distribution function in its flnal form does not carry 
any color degee of freedom. 

We will now set up the relevant equations with the color degree of freedom, following several 
works of Heinz[29]. In non-abelian theory the color charge is a continuously var\'ing function of 
time. The precession of the color charge(Q“) obeys Wong’s equation[30]- 


dr 


( 2 . 12 ) 


which supplements the Lorentz force equation 


dp^ 
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(2.13) 


Here is the gauge potential, and is the structure constant of the gauge group. In order 
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to include the color charge m the phase-space, we consider an extended one particle phase space 
of dimension d = 6 -f (/V- — 1) m SU(N) The extended phase space is taken to be the direct 
sum if® © G, where G is the( compact) space corresopndmg to the guen gauge group In short, 
in addition to the usual 6 dimensional phase space of coordinates and momenta we now have 
another eight coordinates corresponding to the eight color charges in SU(3). In this extended 
phase space a typical transport equation reads as[31]' 


'y,d^ + -f- 


f{x,p,Q) = C{x,p,Q) + S{x,p,Q) 


( 2 . 14 ) 


Here f{x,p,Q) is the single particle distribution function in the extended phase space The 
above equation is gauge and Lorentz in\ariant[31]. The first term in the left hand side of 
equation (2.14) corresponds to the usual convective fiow, the second term is the non-.A.belian 
version of the Lorentz force term and the last term corresponds to the precession of the charge 
as described by Wong's equation. Note that, in general we have to write separate equations for 
quarks, antiquarks and gluons since they belong to different representations of the gauge group. 
For anti-quarks the distribution function f{x^p,Q) obeys a simillar equation, with replaced 
by —Q°' {i.e. the second term in the above equation changes sign). These equations are closed 
with the Yang-Mills equation, given by 

{Df,F‘^'')^{x) = f{x) = g J p''Q'^[fg(x,p,Q) - fg{x,p,Q) + fg{x.p,Q)]dPdQ (2.15) 
The integration measures in the above equation are 

dP=2e(p„)%=-m“)^ = -^ (2.16) 

(27r)'^ [ZTTypo 


dQ = SiQ’^Q^ - q^idabcQ^^Q’^Q'^ - f)d^Q- 


(2.17) 
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Using these distribution functions we obtain the local number density and local energy density 
of the quark-gluon plasma by the equations 

n{r) = J dPdQ{p^u^){fg{x,p,Q) + fg{x,p,Q) + fgix,p,Q)) (2.18) 

€(r) = y dPdQ{p^u^)^{fg{x,p,Q) + fq{x,p,Q) + fg{x,p,Q)) (2.19) 

In the next chapter we shall consider the simplest of the non-.A.belian group viz SU(2) to highlight 
the dynamical role of the color degree of freedom. The realistic guage group SU(3) applicable for 
QGP will be studied in detail in chapter 4 Both the chapters contain rather ad hoc adaptation 
of the Schwinger formola for q^ and gg production that will be reconsidered m chapter six. A 
full study of the transport equation with the new source term is beyond the scope of this thesis 
That will be reported elsewhere. 



Chapter 3 


Solution of transport equations and 
evolution of qq plasma 

3.1 Transport equations for qq plasma 

In this chapter we discuss the production and the evolution of quark-antiquark plasma in URHIC 
by employing the Flux-tube Model[18, 19], m particular to obtain a time devlopment of the quark 
and antiquark distribution function. As mentioned earlier, in this model the two nuclei that 
undergo a central collision at ultra high energies are lorentz contracted as thin plates[21, 22]. 
These two lorentz contracted nuclei pass through each other and, in the process, acquire a non- 
zero color charge {{Q) = 0, {Q'^) ^ 0), by exchanging soft gluons. So one may figuratively call 
such nuclei after collision as color capacitor plates between which a strong chromoelectric field 
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is created. This field creates and gluon pairs via the Schwinger mechanism which enforces the 
instability of the vacuum in the presence of an external electric field The partons so produced 
will collide with each other and also get accelerated by the parent background field. The mutual 
collisions drive the system towards equilibrium with suitable modulation from the background 
acceleration. In all the earlier studies (e g [19, 25]) a set of transport equations are solved within 
this model in the abelian approximation, i.e. the model was valid for an electromegnetic plasma. 
We now attempt to study the evolution of a non-Abelian plasma by solving transport equations 
which are written in an extended phase space of coordinate, momenta and color. 

As mentioned eariler, for simplicity and also to get a feeling for how crucial is the non-abelian 
effect in the study of quark-gluon plasma can be, we shall study the gauge group SU(2). and 
consider only a qg plasma. We will solve the transport equations in SU(3) in the next chapter 
with gluons taken into consideration. For SU(2) charges that we are interested in. the non- 
Abelian extension of Eqn.(2.14) reads 


-b Q^F^yd; + 


f{x,p,Q) = C{x,p,Q) + S{x.p,Q). 


(3.1) 


Note that we do not have to write a separate equation for anti-quarks since they belong to the 
same color representation as quarks. The quarks considered in this chapter does not possess any 
other quantum number such as the electric charge. In reality the situation is of course different. 
Indeed, the antipodal points on the sphere (corresponding to the color part of the phase-space 
with a fixed charge) represent the particle and anti-particle. Wong’s equation (Eqn 2.12) may 
be recast in the more convenient form 
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dQ 

dr 


= u^Q X 


(3.2) 


where the arrows now denote the direction m the color space. In the same notation, the transport 
equation obtains the form 


f{x,p,Q) = C{x,p,Q) + Six.p.Q). (3.3) 

For the model at hand, let us take the ‘plates’ to be moving along the z-direction. Also we 
restrict to only that field configuration for which there exists a gauge choice such that the 
gauge potentials commute with each other every where. Keeping this in mind, we make the 
gauge choice where only the components .4^“ = (A°^, are non-vanishing and the other 
components are zero. With this choice, the resulting electric field points in the ‘3’ direction in 
the color space. We next require a boost invariant description [2] of the distribution functions 
as well as the other physical quantities which may be determined thereof, even as the system is 
evolving. We now proceed to impose the Lorentz gauge condition, which has to be done such 
that the chromoelectnc field depends on the boost-invariant(along the axis of collision) quantity 
r — (^2 _ 2 : 2 )i /2 Observing that the system is effectively (1 + 1) dimensional, in the t-z plane, 
we write 


+ Q ■ F,,p^d; + p^Q X A, ■ 


dQ 


A'^ = t^^d,G{r) (3.4) 

where the indices p, v take values 0, 3 and G is a Lorentz scalar function depending only on r. 
The above choice automatically implements the Lorentz gauge condition dp,A^ = 0. Clearly, the 
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chromoelectric field E is dependent only on r, and is given by 

G(t) (3.5) 

(It Will be seen below that all the derivation in this chapter is obtained in terms of E(r), not 
G(r). However this relation (Eqn. 3.5) is important (see chapter 6) when we use new source 
term (from r dependent electric field) for particle production.) 

The Wong equation (Eqn 3.2) guarantees the conservation of the magnitude of the vector charge 
Q, which may now be held fixed. Being the analogue of the Larmor equation for a charged 
particle in an external magnetic field, it also conserves the component of the charge that is 
parallel to the background chromoelectric field. It is, therefore, convenient to resolve the SU(2) 
charge in the polar coordinates. Writing 


E{r) = 


dr ! \ ^ 


Qx = Qsindcoscf)] Q 2 = QsinOsincj}-, Qz = Qcosd , 


(3.6) 


it is straight forward to verify that 


A.\ 

dQ;, s<i, 


(3.7) 


The last equation will lead to considerable simpification in solving the transport equation. 


Let us now consider the collision term. The color plate model can be looked upon as an effective 
version of a more complete microscopic theory describing the interactions of quarks and gluons. 
Hence, we employ the simple relaxation time hypothesis, with a phenomenological parameter 
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Tc, the relaxation time, to be obtained from microscopic computations. It then follows that 


C = 


- feq) 

Tc 


(3 8) 


where /e, is the eqilibrium distribution function, with local (space time dependent) values of the 
thermodynamic quantities. While it has been customary to take /g, to be one of an ideal gas, 
with some support from the parton cascade model (pcm) analysis of Geiger and Kapusta[17], 
recent lattice studies suggest that the quark-gluon phase is possibly a non-perturbative phase, 
being ideal only at T ~ STg. Since we are not dealing with a true system here, we shall 
conveniently take /e, to be that of an ideal Fermi gas with a local temperature, evolving as 
a function of the proper time r Apart from space time dependence this distribution function 
carries color in its definition [32]. We thus have 


/e,= 


2 

exp((p^ — A'")u^/T(r)) + 1 ’ 


where is the flow velocity. 


(3.9) 


u’^ = {coshri^0,0,sinhT]). (3.10) 

which is written in terms of the (space-time) rapidity tanhrj = z/t. 

It can be verified that with the above choice of the potential 


= 0 , 


(3.11) 
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and the above equation boils down to the familiar form 

2 

exp(p^)u^/r(T)) + 1 

Now we demand the boost invariance following Bjorken’s picture according to which the lon- 
gitudinal boosts are the symmetry operations on the single particle distribution. The boost 
invariant parameters on which / can depend are, apart from the charge coordinates, 

T = ( = {ri - y), pt = {pi - pIY'^ (3.13) 

where y = tanh“^(p;/po) is the momentum rapidity. It is convenient to write separate equations 
for quarks and anti-quarks, a la the Abelian case. If we therefore identify the quark states 
with the points on the upper hemisphere of the color sphere, and the antiquarks with the lower 
hemisphere, by a trivial relabeling, we may write two equations. 


(3.12) 


dr 



+ 


g cos 0E{r) 
Pt cosh ^ 


d 


+ g—G{r) tanh^ 
dr 


d 

d(j> 


f{T,(,pt,9.6) 


+L = 

Tc Tc Pt cosh ^ 


(.3.14) 


d f tanh ^ 

5r \ r 


g cos $E{t)\ 
Pt cosh ^ / 


d 


d 

+ g—G{T) tanh^ 
dr 


d 

dcf> 




Tc Tc Pt cosh if 


(3.15) 


with the first of them for the quarks and the next for the anti-quarks. The angle variable 9 
varies from 0 to 7r/2 in both the equations. 
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Finally, S is the (non perturbative) Schwinger’s expression for pair production and is given by 




gE cos d 
Stt^ 


In 


— exp 


27rp? 


gE cos 6 


(3)^/^exp(-Q<f^) 


(3.16) 


where we have inserted the Gaussian dependence on ^ by hand. Note that the non-Abehan 
nature features dominantly ma the cos9 term This occurence plays an important role in the 
evolution of the system. 


After having set up the relevant equations (not all yet since energy-momentum conservation is 
to be imposed), we observe that the above differential equation possesses the formal solution 


= / drexp( 
Jo 



s(r^e^p^g) ^ /e.(T^e^Pt) 

Pt cosh E '^c 


where <^(r/) is given by 


(3.17) 


(f' = smh ^ 


— sinh + 


g cos 6 

PtT' 



dr"E{T") 


(3.18) 


It is thus clear that / does not depend on d, so the ^ term contributes nothing to the transport 
equation. However the 9 dependence is still involved through out the formalism which shows 
the non-Abehan effects. The corresponding distribution / for antiquark can be written by just 
changing gto—gmf as can be checked out easily 


The transport equation has to be supplemented with the conservation constraint 
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d,TZt + = 0 , 


(3 19) 


since it is the field energy that is being pumped in order to produce the qq pairs More explicitly, 


d,TZ, = -fJZ s -d,TZ, {3.201 

with 

TZ,= jpVU + DdTdP-Q (3.21) 

where the measures dO.Q = sm 9 ddd(f> , dP = Ja — I P^Q^if ~ f)dTdClQ. 

Here we have taken the value of the degeneracy factor 7=2 corresponding to two flavours. 
Since energy and momentum are conserved in each collision, the moment of the sum of the 
colli.sion term vanishes. 


J p"(C)drd.Qg = 0 (3.22) 

Now taking the first moments of the Boltzmann equation and integrating over the color degrees 
of freedom for / and / and making use of the conservation of energy and momentum, we obtain 


djr + gEir) / dVdO.Qp^ 


3 ( 


+ 2 / dTdnciS = 0 


(3.23) 


where 
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T;'' = (iiag(£V2, E'^12, EV2, -E'^ 12) (3 24) 

IS the energy momentum tensor for the field We may solve for the electric field by employing 
the same procedure as in[25], and by employing the symmetry /(r,if.px,0) = /(t, — 0). 
We thus obtain the equation governing the decay of the source field to be 


dE{T) 

dr 


lo (f6'sin6'[/ - /] + Ya|£'(r)|^/^ = 0 


(3.25) 


where a = aC(5/2) exp(0.25/o;), a = c{gj2)^^‘^j^^ and c = Finally, C(5/2) = 1.342 is the 
Reimann zeta function. 


Equations (3.25) and (3 17) are as yet underdetermined since the local temperature T(r) is free. 
In order to fi.x the form of T(r), we appeal to the relaxation time approach that we are employing, 
and assume that the the particle energy density differs negligibly from the equilibrium energy 
density. We may then relate, by an ansatz, the proper energy density which is defined by 


£(r) = J c(ro(ng(p»u„)2(/ + /) (3.26) 

to the temperature by its equilibrium value, whence, 

(3 27) 

It may be mentioned that the weaker condition of energy-momentum conservation that we have 
employed here , in fact, satisfies the Yang-Mills equations ^ as well. 



28 


Finally, we pause brief!}’ to discuss the effects of hard thermal loops that have been emphasized 
recently [33]. They have been derived from the classical transport equation as well by Kelly 
et a/[34, 35] The latter derivation, which is important for us here, is based on an analysis of 
the Vlasov equation which properly describes the expansion of an already equilibrated gas, in 
a back-ground field but with with no source or sink. It is. therefore, necessary to determine 
whether the more complete transport equation such as the one that we have here will have any 
effect on the results of Ref. [34, 35] Conversely, it is also necessary to determine how the hard 
thermal loops will affect the purely classical results that will be obtained here. Let us recall that 
the derivation of Kelly et al[34, 35] consists of a systematic expansion of the distribution function 
as well as the background field term in powers of the coupling constant. The zeroth order term 
is merely the Fermi-Dirac term for the quarks. Since we are solving here the approach to the 
equlibrium, there is necessarily a dependence of the distribution functions and other physical 
quantities on g. It remains to disentangle the contribution coming from the hard thermal loops 
Indeed, observe that the two extra terms that we have at hand here are of higher order in g. 
First of all, the source term is non-analytic and does not even admit an expansion in powers 
of g The other collision term is easily seen to be of order g'^ or higher. We thus conclude that 
the conclusions of Ref. [34, 35] remain unaffected, and we may take over those results in toto, 
supplementing our classical results 


3.2 Computational Procedure 


We have adopted here a double-self consistent method to determine (/, e, ...) , following the work 
of Banerjee et al[25] in the Abelian case. The procedure follows the scheme {T{r)triai, F^('r)trtai} — t 
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{/, / E{t)} {/, /} T{t) -> ... by repeated use of equations (3.17), (3.25), (3.17), (3.27). 

The iteration terminates as soon as a convergence is established in the solutions for £(r).r(r) 
All the desired quantities are thereby consistently determined. The same procedure will be 
employed subsequently. 


3.3 Results and discussions 

Before we proceed to present and discuss the results, a few comments about the choice of 
the value of the parameters in the computation We put g = 4 throughout our calculations. 
Since lattice computation results predict a phase transition from the baryonic phase to the 
quark-gluon phase at densities ~ 5 — lOGeV/ it has also been customary to take an initial 
energy density in the same range m the Flux tube model. However, greater care needs to be 
taken before the initial energy densities in URHIC are chosen. Indeed, a fairly reliable estimate 
of the time required for the partons to be produced in central collisions after the two nuclei have 
suffered the maximum overlap is ~ 0.05 —0.1 /m; a simple dimensional analysis leads to a value 
of the initial energy density to be {1/2)Eq « 500 — 1000 GeVj fm? ^ the precise value depending 
on the the magnitude of g in Since in the color plate model, we set the zero of 

time not at maximum overlap, but at the instant when the plates have completely crossed each 
other, we shall take the initial energy density e = ZOOGeVj fm^ . We shall also be guided by the 
results of Geiger and Kapusta[17i m our choice of the values of Tc, and take Tc = 0.2/m. to be a 
realistic value[17]. The scale for the hydrodynamic limit will be set by the choice for as well 
as the formation time tq. A typical value for the analysis at hand is 0.001 fm, which we employ 
here. Finally, we study the other extreme ctise, the collisionless limit, by the choice Tc = 5/m. 
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It IS instructive to compare how the realistic regime behaves in relation to the results that we 
obtain for the two extreme limits. 

We have studied the decay of the chromoelectric field, the evolution of the particle energy and 
number densities, the evolution of the distribution function and its approach to the equlibrium 
state We have compared them with the corresponding Abelian results. In addition, we ha\'e 
also evaluated a quintessentially non-Abehan quantity, viz, the expectation value of the angle 
that the quark charge makes with the direction of the chromoelectric field. In the current model 
the chromoelectric field not only decays to produce the qq pairs, but is itself built up by the 
receding nuclei which act as color plates. We have, therefore, also studied the particle energy 
per unit transverse area as a function of ordinary time to highlight this feature. 

3.3.1 Discussion of the results 

We shall present the results of our analysis in the three regimes corresponding to the hydro- 
dynamic, the realistic, and the collisionless cases. Comparison between the .Abelian and the 
non-.A,belian systems will be taken up subsequently. 

Consider the hydrodynamic limit first. The physical quantity of utmost importance is the 
chromoelectric field, whose dynamics is most readily determined in the Flux tube model at 
hand Indeed, apart from determining the production and the acceleration of the partons m the 
Flux tube model at hand, it has an additional role which has been emphasized by Svetitsky[36]: 
the charms which are produced in the pre-equilibrium stage would not only interact with the gas 
of light quarks and gluons, before either forming a J/i/; or an open charm mesonic state, but will 
also be influenced by the background field. Thus the study of the evolution of the mean electric 
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field in the pre-equilibnum stage acquires an added importance. It is a merit of the Flux tube 
model that vve can readily determine the evolution of the background electric field. Note that in 
contrast, no such information can be extracted in the more microscopic models such as parton 
cascade model (PCM)[17] and heavy-ion jet interaction generator (HIJING)[37h In Fig 3 1 the 
decay of the field is shown Recall that = .001/m As can be seen from Fig 3 1. the electric 
field has hardly decayed at all, with a percentage decay less than 2 %, even at r = 1 5/m. 
The corresponding particle energy density, shown in Fig. 3.2 is also negligibly small, with the 



Figure 3.1; Decay of the chromoelectric field as a function of proper time (in units of fermi), for 

Tc — .001/m The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

ratio with respect to the initial energy density being ~ lO""*. The Abelian situation is an 
order of magnitude better, but is still hopelessly small. And indeed, the number density is also 
very small, and as shown in Fig. 3.3 gets saturated at ~ .Oof fm^ and the temperature also 
stabilizes to a value ~ 40Mey, which is much less than the temperatures required. Clearly, it 
is very unlikely that the plasma does not go through a prequilibrium phase. 


The behaviour of the system for = 0.2 and 5/m is in sharp contrast to the case discussed 
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T (fm) 

Figure 3.2: The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for = 001/m. The solid line refers to the non-Abelian 
case, and the broken line to the Abelian case. 

above There is a significant decay of the electric field, which is shown in Figs. 3.4 and 3.5 
for the two respective values of r^. The electric field decays by about 15% of the original 
value, at r = 1 5/m. The corresponding energy densities are also not very different, and 
yield ~ 10 — 15% of the initial field energy density (see Figs. 3.6 and 3.7). The corresponding 
temperatures are ~ 300Me%, which is quite realistic. Note that the parton cascade model[17] 
predicts a value about 300 MeV at 2.4/m. The real difference between the collisionless case 
and the ’’realistic” case is in the number density. Indeed as Figs. 3.8 and 3.9 show, the plasma 
produced for = 0.2/m has a number density ~ 15 — 20//m^, the corresponding number for 
Tc = 5/m is three to four times smaller. In other words, the system is more dense in the former 
case than the latter. Clearly, this distinction should show up in signatures such as dilepton 
production which are sensitive to both the number as well as the energy density. 

It thus follows that the flux tube model, with the incorpration of the color degrees of free- 
dom, definitely rules out instantaneous equilibration as envisaged in the approach of Baym[18]; 
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t(fm) 

Figure 3.3: The particle number density as a function of proper time (in units of fermi), for 

Tr = .001/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

further, it also distinguishes the collisionless limit from realistic values of equilibration time. 

Finally, it should be emphasized that the flux tube model as employed here not only pumps 
in the held energy to particle production, but also contributes to the fleld energy by virtue 
of the recession of the color charged nuclei from each other. In fact, a crude estimate of the 
time required to convert all the plate energy to the fleld energy turns out to be ~ 5/m for 
200Gey/nucleon In any case, it is therefore misleading to interpret the ratios we have shown in 
Figs. 3. 2, 3. 6 and 3.7 as the fraction of the total fleld energy that has gone into the particles. To 
emphasize this, we have evaluated the dependence of the field and particle energies per unit area 
by integrating over the contribution along the longitudinal direction. The results are shown in 
Fig. 3.10, from which it is clear that there is a lot more energy m the field than in the particles. 
For the same reason, it is also misleading to interpret the field energy density at r = 0 as the 
counterpart of the energy densities employed in lattice computations to study the transition 
from the hadronc to the QGP phase. 
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Figure 3.4: Decay of the chromoelectric field as a function of proper time (in units of fermi), 

for Tc = 2/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian 
case. 

3.3.2 Comparison with Abelian results 

It IS clear from Figs. 3. 1-3.9 that the Abelian and the non-Abelian results bear little resemblence, 
belying the expectation that the ’’Abelian dominance” holds in this case The difference is most 
dramatically highlighted at = 0 2/m. Whereas almost all the initial field has decayed in the 
Abelian case, only 30% has done so in the non-Abelian case. Accordingly, The particle energy 
density is larger by a factor of ~ 5 , the number density by a factor of ~ 4, and yields an 
abnormally large value of~ 800 A/eV" for the temperature. Earlier Abelian calculations[19, 25, 
28] yielded reasonable values for the temperature because of unrealistic values for the initial 
field energy It is a universal feature that the non-Abelian plasma is rarer and cooler than its 
Abelian counterpart. In the hydrodynamic limit, the Abelian ananlysis yields a temperature 
~ 200MeV, in contrast to ~ 40MeV in the non-Abelian case. Thus the incorporation of the 
color degree of freedom rules out instantaneous hydrodynamic evolution. Of course, a colorless 
plasma does not have reasonable temperature for any other r.-. 
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Figure 3.5: Decay of the chromoelectric field as a functioa of proper time (in units of fermi), 

for Tc = 5/m. The solid line refers to the non- Abelian case, and the broken line to the Abelian 
case. 

Yet another interesting aspect that emerges from our studies is in the close interplay between 
the value of rc and the color degree of freedom While the field decays faster in the non-.A.belian 
case in the hydrodynamic limit, the trend reverses at = 0.2/m and gets restored in the 
collisionless limit. In contrast, the particle energy density is smaller for a colored plasma in 
hydrodynamic limit. Further, it continues to be so. at = 0.2/m but becomes larger than the 
Abelian case in the collisionless limit. These non-trivial manifestations of the color charge may 
be expected to have an important bearing on the other bulk properties of the plasma. 

In order to gain some insight into the above features, we have also evaluated the expectation 
value of 0 =< cos^O > which yields the rms value that the particle charge makes with the field 
direction in the color space. Note that the above quantity is gauge invariant, and hence physical. 
Fig. 3.11 shows 0(rj for = 0.2/m. It may be seen that the value saturates around 0.25, 
corresponding to 0 ~ r/3. The corresponding value in the non- Abelian case is strictly zero. 
The effect of the background field is thus reduced, leading to dominance of the collision 
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X (fm) 

Figure 3 6: The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for = 2/m. The solid line refers to the non-Abelian case, 
and the broken line to the Abelian case. 



X (fra) 

Figure 3.7; The particle energy density scaled w r.t the initial field energy density as a function 
of proper time (in units of fermi), for — 5/m. The solid line refers to the non-Abelian case, 
and the broken line to the Abelian case. 
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Figure 3.8: The particle number density as a function of proper time (in units of fermi), for 

Tc = 2/m. The solid line refers to the non- Abelian case, and the broken line to the Abelian 
case. 



t(fm) 

Figure 3.9: The particle number density as a function of proper time (in units of fermi), for 
Tc = 5/m. The solid line refers to the non-Abelian case, and the broken line to the Abelian case. 
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term in the expansion of the particles Indeed, it shows up most clearly in the approach to 
the equilibrium, where one now expects that the plasma equilibrates the fastest in the direction 
normal to that of the field (in the color space) This is corroborated as may be seen in Fig. 
3.12 



t(fm) 

Figure 3.12: f j feq as a function of proper time at Pf = 200MeV', (^ = 0 and = 0.2/m for 
three different angles corresponding to cos9 = 0 (solid line), cos6 — .25 (dash line just below the 
solid line), and cosd = 1 (the other dash line). Note that the equlibration is fastest aX 6 = ii 12. 


3.4 CONCLUSION 

To conclude, we have studied the production and the equilibration of a genuinely non-Abelian 
plasma with the color degree of freedom incorporated in both the source and the background 
field term m the transport equation. In the SU{2) gauge theory that we have considered here, 
the distribution function is defined in the extended phase space. We find that this approach 
recaptures in an elegant manner many of the findings of the more microscopic parton cascade 
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model. It has the further advantage that it indeed exhibits the color degree of freedom manifestly, 
and allows us to compute various gauge invariant quantities Significantly, we also find that the 
Abelian approximation, employed hitherto in most studies of equilibration of QGP is rather 
too drastic to be used for any quantitative analysis and comparison with the experiments The 
study also almost rules out instantaneous equilibration, and also strongly suggests that the 
collisionless limit may also not be the favoured in URHIC. 

To be sure, we have not made any comparison with the experimental findings here, for the simple 
reason that we are as yet dealing with a simpler gauge group , and ignored the gluonic component 
altogether. The indications from the present study are unmistakable, though. Indeed, the 
particle production is enhanced because of increase in the phase space available, and for the 
same reason the plasma will be cooler than the Abelian counterpart. The energy goes to the 
color degree of freedom, and does not simply heat the system as it would happen m a colorless 
plasma. If we consider the realistic 517(3) case, this feature will get further accentuated; For 
the same initial configuration, we may expect a rarer and a cooler plasma. Of course, there 
are other features which are intrinsic to 517(3): there is yet another Casimir invariant, and the 
gluon term will also have to be incorporated. They will be considered in the next chapter. 



Chapter 4 


Solution of transport equations and 
evolution of quark-gluon plasma 


4.1 Transport equations for quark-gluon plasma 

In the last chapter we have studied the evolution of a non-abelian qg plasma within the color flux- 
tube model where it was shown that the non-abelian features have a major role in the evolution 
of the system in a manner that was not captured in any of the earlier studies [18. 19, 25]. It was 
pointed out that it was unlikely that the system would equilibrate instantaneously. However, 
for simplicity, we had ignored the gluonic component and considered the simpler gauge group 
SU(2). In this chapter, we remedy both the drawbacks; we consider the gauge group SU(3) 
as is appropriate for a real plasma We also include the gluons, which have been completely 
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ignored so far We pay attention to the detailed evolution, and its prediction for different bulk 
properties of the plasma in this chapter. Signatures such as dilepton production within this 
model will be studied in the next chapter. 

There are models such as HIJING[37] and PCM[38] which study the pre-equlibrium dynamics 
via perturbative QCD (pQCD). The production and the evolution of hard and semi-hard partons 
are studied by a master rate equation in HIJING and by a transport equation in PCM. While 
they may be reliable in the study of hard partons, these perturbative approaches are admittedly 
insufficient to study the dynamics of the soft partons[37], in particular, their production. Being 
perturbative, they do not incorporate any of the non-perturbative aspects such as the formation 
of the strings and their break up, which is studied even in pp collisions. The addition of soft 
pajtons to hard and semi- hard partons changes the bulk properties of plasma, such as temper- 
ature and energy density. Xu et a/.[10] have observed that this addition leads to an enhanced 
suppression of J/tP, which can be understood to be a consequence of increased number density 
and a lowered temperature of QGP. In our study, we employ the Schwinger mechanism[23] for 
particle production, which is quintessentially non-perturbative. Even if one were to employ a 
perturbative version, as a time dependent electric field would require[39], which we shall propose 
in chapter six, it may be noted that an initial electric field as a classical saddle point owes its 
existence to non-perturbative processes, viz, the soft gluon exchanges that take place between 
the two nuclei. On the whole, one may expect that the pQCD based studies[37, 38] will be 
useful in the hard regime, and that the flux tube like models will be required to study the soft 
regime which will become increasingly prominent as the system expands and more and more 
secondaries are produced[40]. Indeed, it should be possible to develop a unified approach to 
study the hard and the soft components, say a la the approach of Eskola and Gyulassy[41] who 
have included the minijets in their so called chromoviscous hydrodynamics which is again based 
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on the flux tube model. This needs a separate study. 

For our purpose here, the transport equations we solve in the SU(3) gauge group to describe 
the evolution of quarks and gluons (see chapter 2 for a detail description) are: 


[p,d^ + 


f{x,p,Q) = C{x,p,Q) + S{x.p.Q). 


(4.1) 


Here f{x,p, Q) is the single particle distribution function in the extended phase space of dimen- 
sion 14 which include eight color coordinates in SU(3) gauge group in addition to the usual six 
dimensional phase space of coordinates and momenta. The first term in the left hand side of 
the above equation corresponds to the usual convective flow, the second term is the non-abelian 
version of the Lorentz force term and the last term corresponds to the precession of the charge as 
described by Wong’s equation. S and C on the right hand side of the above equation correspond 
to the source and collision terms respectively. 

Note that we have to write separate equations for quarks, antiquarks and gluons since they 
belong to different representations of 517(3). 

The term gluonic source merits some explanation here. The claissical background field that we 
consider here has, in contrast to the Maxwell field, self interaction. We are interested in the 
stability of the gluonic vacuum ( which is the analogue of the photon vacuum in eletrodynam- 
ics) against the fluctuations in the classical background field. An adaptation of the Schwinger 
mechanism m QED shows that the fluctuations can indeed produce the gluons, i.e. the off-shell 
classical field can spontaneously produce the on-shell radiative gluonic field (see Eqn 4.7). There 
IS, therefore, no ambiguity or double counting in this process. The source term yields asymptotic 
gluonic states, and further interaction between the gluons is treated separately by a collision 
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term. 

It is very difficult to solve the transport equation written above, in general. And we have a set 
of three coupled equations here As was done in the earlier chapter we make a few assumptions. 
First of all, we admit only those potentials which can be brought to a form where the only 
surviving components are = (A®'', This choice restricts to be “Maxwell” like, 

also pointing in the “1” direction m the color space. This restriction is not entirely arbitrary as 
it is known[42] that the non-Maxwellian configurations - where the charges, the gauge potential 
and the fields do not lie in the same direction in the group space - do not produce particle 
pairs, in general. Secondly, we require a boost invariant description of the physical quantities [2]. 
Accordingly, we demand that the distribution functions also be functions of only boost invariant 
quantities. Finally, we work within the Lorentz gauge which is implemented elegantly by the 
choice A^“ = e^‘'duG°-{r), n, u = 0,3, with all the other components zero, r = {F - is 

the boost invariant proper time. 

Now we fix the magnitude of the vector charge which corresponds to the first Casimir 
invariant of SU(3). It is simply the value of the coupling constant. There is also another 
Casimir invariant, viz. d°-^^Q°-Q^QY which is also conserved as the QGP evolves. There is, 
however, no way of fixing its value and the experiments presumably impose no restriction on its 
allowed values. In fact, the same holds for lattice analyses as well. For this reason, we do not 
take cognizance of this invariant, and we conveniently resolve the SU(3) charges in the polar 
coordinates: Qi = Q ni=\ sindkCosO^, for i 7^ 8, and Qa Q nl=i sindk- 

As is studied by Heiselberg and Wang[4.3] a collision term can indeed be obtained from pQCD. 
Apart from making the equation hopelessly non-linear, this approach would be good only for 
the hard components, which are not of interest to us here. As aptly pointed by Hung and 
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Shuryak[44] recently, a microscopic description of collisions gets increasingly cumbersome and 
also unnecessary as more and more secondaries are produced. On the other hand it is admittedly 
true that there is no way to reliably obtain a collision term in the non-perturbative regime. So 
we shall employ a relaxation time approach here, where the relaxation time will have to be 
fixed phenomenologically, and in all probability, a posteriori. In reference[43] it is proposed 
that Tc can be local and have a (weak) dependence on r. We take to be constant here, and 
further refinements may be incorporated after we have a better understanding of the transport 
phenomena both experimentally and theoretically. 


Within the relaxation approach, the collision term is written as 

^ - f’") 2) 

Here is the distribution function, in local equilibrium. Note that the locality can extend to 
the color space as well apart from space-time. Taking it to be an ideal gas, for simplicity, we 
write, 

2 

^ exp((p^ - Q-A^^<^)uJT{t)) ± 1 

where the -l-(-) sign is to be taken for fermions (bosons). Note that we allow a common equi- 
librium temperature for quarks and gluons. Here = {coshr],0,Q, sinhrj) is the flow velocity 
and rj is the space time rapidity given by tanhrj = zjt. With our choice of potentials, it follows 
that = 0, so that the equilibrium distribution function can be written as: 


req 


exp((p'^u^)/T(r)) ± 1 


(4.4) 


Note also that the temperature depends only on the proper time, in accordance with the Bjorken 
picture. Demanding the same of the distribution functions as well, we require that the longit- 
udinal boosts be symmetry operations on the single particle distribution. As mentioned earlier 
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the boost-invariant parameters are, apart from the color charges, 


r = {f~ ^2)1/2, ^ = {rj- y), pt = {pi - pf) 


1/2 


(4.5) 


where y = tanh ^{pi/po) is the momentum rapidity. 


The above set of invariant variables also serve to express the source terms for qq and gluon 
pairs which will be obtained by the Schwinger mechanism The expression for qq production 
(obtained from constant electric field) is given by[45, 46] 

2irpl 


u i c a \ 9E cos 9i 
^q{TA^?tJi) = 5-5 — In 


1 — exp 


-Cl, 


)^^^ exp(— Qi^^) (4.6) 


TT 


gE cos di 

and for gluon pair production the corresponding term in SU(3) is given by tl j relatively stronger 
term [31, 47] 

S,{T,i,P,M = (3/2)5, (T,f,p„9,). (4.7) 


We take g — A throughout our calculation. 

Now consider the third term Q), in the transport equation (4.1). Due to 

the restriction of the gauge potentials to the form = (A°\A^^), there is an additional 
simplification in the transport equation. In order to see that, note that equation (4.1) hats the 
formal solution 


where ^'{rf) is given by 


dr' exp{- -) 


SqAr\i',Pi.Q) , fzi'^',e,PuQ) 

Pt cosh Tc 


(4.8) 


r' ptT' Jt' ' 


e' = sinh“' 


(4.9) 
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Clearly, f^'^’^A^Q^-^f{x,p,Q) = 0, as S, /«, and only depend on Qi 

With this final simplification, we get a set of three equations, one each for quark, antiquark and 
gluon respectively. They are, explicitly. 


d 

dr 



g cos9iB{t)\ d 
Pt cosh^ J 
f f®'’ 

I -'g'g _ ■'g.g 

Tc Tc 


f <1<9 ( ‘ ^ i Pi ^ ^l) 
Pt cosh ^ 


(4.10) 


for quarks and gluons, and 


d 

dr 



geos 6 iE{t) 
Pt cosh ( 


d 




^ii ^ IL + MhHiMA 

Tc Tc Pt cosh ^ 


(4.11) 


for antiquarks. 

In the process where the field and the particles are present the conservation of energy-momentum 
is expressed in the form: 

d,TZ, + = o. ( 4 . 12 ) 

where 

TZi = / P>''(2/, + 2/, + /,)drdnr, (4.13) 

and 

r;" = diag(£:V2, E^I2, E^/2, -E^/2). (4.14) 

Here dV = d^p/(27r)^po = Ptdptd^ t {2Tr)‘^ , and dfij is the angular integral measure in the color 
space. The factor 2 m the equation (4.13) is for two flavors of quarks. Now since energy and 
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momentum are conserved in each collision, we have: 

J p'^CdVdO.T = 0. (4.15) 

Taking the first moment of the Boltzmann equation, integrating over the color degrees of freedom 
for /,, fq and fg and making use of equations (4.12) and (4 15), we obtain from equation (4.1) 


+ gE{T) j ^4 J dTd^ijp’'Sq + 2J dVd^rp^Sg = 0 (4 le) 

In the above equation the factor 4 in third term arises because we have two separate transport 
equations for quark and antiquark, each coming with two flavors, and the factor 2 in fourth term 
is due to gg pair production, although there is only one transport equation for gluon. Putting u 
= 0 and 3 in (4 16) we get two equations, which then yield the following equation for the decay 
of the electric field: 

rdptp^, Tdesinhe rd^i[2/,-2/, + /,] + (7rV6)alE(r)p/2 = 0. (4.17) 

CLT Ztt Jo Jo Jo 

Here d = aC(5/2) exp(0.25/a:), a = c{gl2)^^‘^j^^ and c = Finally, C(5/2) = 1.342 is the 
Riemann zeta function. 

To solve this equation we fix the form of r(r), by demanding that the particle energy density 
differ negligibly from the equilibrium energy density, in each collision. We then relate the proper 
energy density, which is defined by 

e(r) = I drdnAp^Ug.y{2fq + 2fq + /,), (4.18) 


to the temperature by its equilibrium value, whence. 
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= ( 4 . 19 ) 

/I 

We solve the equation (4 17) numerically to study the evolution of quark-gluon plasma. 

4.2 Numerical procedure 

The numerical procedure is already discussed in the earlier chapter. We use a double self- 

consistent method. The procedure follows the scheme {T(r)f„a/, —>{/)/) E{r)} — > 

{/, /} — )• T(t) — f .. by repeated use of equations (4.8), (4.17), (4.8), (4.19), which is iterated 
until there is a convergence to the required degree of accuracy. We have in mind the LHC 
energies, and we have taken the initial energy density as cq = 30QGeVf fm? . Compared to 
PCM[38] where the initial particle energy density is taken to be around 1300 Gev/fm^ at 
RHIC, our choice of initial field energy density might appear some what low However, our 
choice IS guided by the estimate that the formation time for a qgp is a fraction of a fermi, as we 
explained in the previous chapter. In any case it would not be very appropriate to compare the 
two initial conditions since the ’’initiality" is only in a limited sense. Indeed, it is the energy 
density at the location of the receding color plates, and being well in the fragmentation region at 
all later times after the collision, what matters is the energy in the central region, corresponding 
to larger and larger values of r as the system evolves. It will be seen that in this region the 
results that our study yields are not unreasonable, if we make a judicious choice for the value of 
Tc. We have studied, in this paper, the results for three different values of For hydrodynamic 
and collisionless limits we have choosen r. = 0.001 fm and 5.0 fm. We have compared the 
results of these limitng value of to a realistic intermediate value = 0.2/m. 
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4.3 RESULTS AND DISCUSSIONS 

The solution of the transport equations following the procedure outlined in the previous section 
allows us to determine directly the temporal evolution of the particle and energy densities, 
temperature and also the rate at which the field energy flows into the particle sector. These 
quantities are of intrinsic interest and are, in principle, amenable to experimental study via 
dilepton production. We present the results below. Also of interest are the broader questions; 
when the equilibration sets in, at what time the flow becomes hydrodynamic - with or without 
viscous flow[2, 41], etc. 

Then there are features peculiar to our model. Since the color plates are receding away from 
each other, energy is continuously pumped into the field, which subsequently decays to produce 
particles. We study the relative rates at which these two dynamical processes proceed in URHIC. 
This process cannot of course proceed indefinitely. This has already been observed by Gyulassy 
and Csernai in their study of the dynamics of fragmentation region [48]. As the plates give up 
their energy, they decelerate. The deceleration sets limits on the times up to which the model 
is valid. Indeed, the assumption of boost invariance breaks down as it is strictly exact only if 
Vpiate = c. It IS safe to assume boost invariance so long as Vpiate > 0.9c[49] A simple estimate 
shows that this condition, for our choice of initial energy density, holds up to ~ 5 — 10/m for 
LHC energies. There would be other competing processes in the fragmentation region, and the 
Bjorken scenario is probably valid upto 3 — 4/m. Keeping this in mind, we have restricted 
ourselves to (proper) times < 1.5/m. Note that the results presented below will be applicable 
for RHIC energies only if r < 1/m. 

Finally, the formulation allows us to study an inherently non-abelian quantity, viz, 



51 


< cos^9 >, where 9 is the angle between the charge and the field in the group space. Note that 
its abelian counterpart = 1. Wherever possible, we have compared our results with the earlier 
SU{2) study and PCM. We also display the distribution function(s) in the color space, which 
none of the other models can yield, be they flux tube based or pQCD based. All the results 
will be shown for the quarks and the gluons separately. We mention here that it is the color 
Flux-tube model within which one obtains the evolution of the mean background chromoelectric 
field, which is absent, either in HIJING or PCM. The evolution of such a background field has 
a greater impact on the acceleration of the partons present in the system, and hence on the 
observed signatures. For example the cc pair, which is produced in early collisions is acted by 
this back ground field along with cc potential, to evolve into a physical J/i/i 

4.3.1 Comparison with SU(2) results 

In the earlier chapter we had argued that the assumption of the so called abelian dominance 
which was made in a large number of transport studies[18, 19, 25] does not receive any justi- 
fication from a proper study of a non-abelian transport equation. The present study reinforces 
the same idea, as it is only to be expected. In particular, note that a determination of quantities 
such as the vacuum current or the polarization current[41] will be particularly suspect in view 
of the fact that the effective charge will only be a fraction of the true charge (see below). There 
now arises the question of the dependence of the results on N if one employs the gauge group 
SU(N) in solving equation (4.1). Although we are not in a position to make a strict comparison 
between SU(2) and SU(3) in this chapter because the SU(2) study had ignored the gluonic 
terms altogether, it is still useful to see what a limited comparison can yield. We shall restrict 
ourselves to the value Tc = 0.2 fm to contrast very briefly the results of SU(2) and SU(3). 
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We shall consider the time dependence of the electric field, the energy.' density and the number 
density, shown in Figs 4. 1-4. 3. First of all, it may be noted that the SU(3) quantities evolve 



T (fm) 

Figure 4.1: Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for Tc = .2/m. The solid line refers to SU(3) case, and the dashed line to SU(2) case. 

much more rapidly than their SU(2) counterparts In fact, while the number density and the 
energy density have attained their maximum value around 1/m, at which time the field also has 
considerably decayed, there is hardly any activity in the SU(2) case even upto 1 5 fm. Indeed, 
at r ~ 1/m, the number density for SU(3) is larger by a factor ~ 10, and the energy, by a 
factor ~ 5, even if we consider only the quark sector. Please note that it is not that the energy 
is merely apportioned between quarks and gluons; the existence of a second channel has in no 
way decreased the flow to the quark sector. 

More significant is the fact that it is the SU(3) flow that shows the required trend towards a 
hydrodynamic flow. The SU(2) coumterparts fail to show any such trend all that way upto 1.5 
fm. We may certainly expect a hydrodynamic flow to occur at much later times, but in all 
likelihood it does not seem to happen at any realistic value of r. We shall discuss the SU(3) 
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T(fm) 

Figure 4.2. The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for = .2/m. The solid line refers to SU(3) case, and the 
dashed line to SU(2) case. 



T(fm) 


Figure 4.3: The particle number density as a function of proper time (in units of fermi), for 
Tc = 2/m. The solid line refers to SU(3) case, and the dashed line to SU(2) case. 





54 


flow in more detail in the next section. On the whole, the plasma is denser for SU(3) due to 
increased volume in phase space. 


4.3.2 Discussion of the results 

We now discuss the results for SU(3) in some detail, for three values of - 0.001, 0.2 and 5 
fm. The first (last) choice corresponds to the instantaneous hydrodynamic (collisionless) case, 
m the time scale set by the initial energy density. The quark and the gluonic contributions will 
be shown separately. It is clear from Figs. 4.4 and 4.5 that the quarks and the gluons show 



X (fin) 

Figure 4.4 The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for = .2/m. The solid line refers to total energy density, 
upper dashed line to quark plus antiquark energy density, and lower dashed line to the gluon 
energy density. 

the same behaviour regarding the energy density and the number density as functions of proper 
time. The quark contribution is larger, partly because we have considered two flavors, although 
the source term[47] favors a larger rate for gluons in the color space. The situation will probably 
get reversed if we include the more correct perturbative source term for quark[39] and gluon 
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Figure 4.5. The particle number density as a function of proper time (in units of fermi), for 
Tf. = 2/m. The solid line refers to total number density, upper dashed line to quark plus 
antiquark number density, and lower dashed line to the gluon number density. 

production[50]. In that case, the single and three gluon production rates are of the same order 
as that of two gluons That the two sectors will continue to show the same trend may be reliably 
assumed. 

The important common feature that Figs. 4.1, 4.4 and 4.5 show is the nature of the evolution 
at times later than 1 fm. The curves suggest an approach to the hydrodynamic flow. Note 
that a hydrodynamic flow would imply that n(r), e(r) and T(r) behave like and r~T3 

respectively. We find that the corresponding exponents are —0.7, —1.23 and -0.31 respectively. 
While the temperature scaling suggests a close approach to the free flow regime, the other two 
exponents show the presence of drag[51], implying that collisions have not completely ceased. 
It may be expected that full hydrodynamic flow will take over around 2/m. 

First of all let us compare the results for quarks and gluons separately, for Tc = 0.2/m before 
presenting the results for different values of r^. In Fig-4.4 we have presented the scaled particle 
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energy densites. In contrast to HIJING and PCM, we obtain lower values for energy and number 
densities for gluons than that of quarks and antiquarks together (Fig-4 5). For a complete study, 
a source term for hard parton production is also required and this term can be obtained from 
minijet production at these collider energies, following say reference[41]. The importance of 
the non-perturbative contribution may be gauged qualitatively by noting that the peak number 
density for quark-antiquark together(for two flavors) in this model is ~ 60 //m^ in contrast to 
the PCM value of ~ 120 j Jtv? (for three massless flavors). The energy density for quark sectors 
in PCM IS also larger by the same factor, implying that the average energy per particle is of the 
same order in both the cases. 

Let us now consider < cos'Q > (where 0 is the angle between color charge and the chromoelectric 
held in color space). This quantity is a good bench mark to characterise the ‘non-abelan ness of 
the system, and we present our results in Fig-4 6. This quantity is gauge invariant and physical. 



x(tm) 

Figure 4 6: < co.s"i9 > as a function of proper time (in fermi) at = .2/rn for quark (solid line) 
and gluon (dashed line) 

It may be seen that this value saturates at ^ m 7r/4, for both quarks and gluons. This value was 
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always unity in an abelian plasma, as 0 = 0 in that case. 

This has a direct impact on the equilibration of the plasma. The equilibration is faster around 
d ~ 7r/2 where the background field effect is zero, and is slower at 0 ~ 0 This is under- 
stood as follows. At larger value of 5, say around 7r/2, there is no acceleration of partons by 
the background electric field. Only collisions are present at this angle, and hence the rate of 
equilibration is faster. On the otherhand when the angle 0 = 0, the acceleration of partons by 
background field retards the equilibration of the partons. At any intermediate values of 0, the 
rate of equilibration lies in between these two extreme values. So this average value, which is 
purely due to the non-abelian nature of the theory, has a major role in the equilibration of the 
plasma. This is more clearly shown in the distribution functions which carry the color degrees 
of freedom explicitly(Fig-4.7 and Fig-4.8). 



Figure 4 7; /a//®’ as a function of proper time at p* = SOOMeV, cf = 0 and = 0.2/m, for 
different values of 0. Solid line corresponds to 0 = 7r/2. 


We conclude the discussion at = 0.2/m by making an interesting observation. .As can be seen 


in Fig-4.6, the x-aiue of < c.o.s‘^9 > fluctutates around its mean value for both quarks and gluons. 
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x(fm) 

Figure 4.8: /g//|« as a function of proper time at pt = SQOMeV, ^ = 0 and Tc = 0.2/m. for 
different values of 9 Solid line corresponds to 9 = it 12. 

with similar fluctuations in other quantities as well - although it is not that pronounced in them. 
While it certainly indicates that the distribution is approaching its equilibrium value asymp- 
totically, the question is whether the fluctuations would persist even as the system hadronises. 
Unfortunately, it is difficult to make any definite assertion at this stage. Indeed, as pointed 
earlier, evolution of the system beyond 3 — 4/m requires a full treatment that goes beyond the 
assumption of the Bjorken scenario in the central region. It is also not known when the hadron- 
ization exactly sets in. However, if one assumes that the fluctuations do indeed persist, it may 
quite well happen that it will manifest as an anisotropy in the parton momentum distribution. 
This has been pointed out in a study of the related ‘color filamentation’ by kIrowczynski[52] 
It is also possible[53], although we do not know in what manner, that these fluctuations in 
the hadronic phase show up as disoriented chiral condensates. It is only a more complete and 
rigorous study that can settle the status of these speculations. 


Now we present the data for the other two values of Tc, namely Tc = 0.001/m and = 5/m, 
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corresponding to hydrodynamic and collisionless limits respectively. These results are then 
compared with = 0.2/m. Consider first the number density n(r). This quantity is seen to be 
sensitive to the values of r^. As can be seen from Fig-4 9, n(r) is larger m the hydrodynamic 



Figure 4.9: The particle number density as a function of proper time (in units of fermi), for 
Tc = 2/m (solid line), for = 0 001/m (upper dashed line) and for = 5/m (lower dashed 
line). 

limit, where the maximum value is around 120 per /m^, at r ~ 1.0/m. As is increased to 
0.2 fm we find a lesser value, the maximum value being around 80//m^. This trend continues 
for collosionless limit where particle density is still less. The increase in number density as 
is decreased is expected in the context of classical theory. This is because, the collision time 
in any classical non-equilibrium theory depends roughly on the inverse of the number density, 
apart from the other factors. The behaviour of number densities in Fig-4.9 for different values 
of Tc in our calculation also reflects the above fact. 

In Fig-4.10 we present the scaled energy densities for different values of Tc- Unlike the number 
density, there is no strong dependence on r,,. This means that the average energy per particle is 
different for different values of r^. VVe find that for r, = 0.001/m, the maximum average energy 
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Figure 4.10; The particle energy density scaled w.r.t the initial field energy density as a function 
of proper time (in units of fermi), for Tc = .2/m (solid line), for = 0 001/m (upper dashed 
line) and for = 5/m (lower dashed line). 

per particle is around 1.0 GeV, whereas it is around 2.0 GeV for Tc = 0.2/m and 4.0 GeV for 
Tc = 5.0/m. Note that these high energy deconfined partons can produce secondary partons, 
such as strange quarks and also can break a fully formed Jj'ip as analysed by short-distance 
QCD In reference[54]. 

Let us consider the electric field. It may be seen from Flg-4.11 that the decay of the field is 
very slow at = 0.001/m, compared to = 0.2 and 5 fm. Only 30 percent of the field 
has decayed in the hydrodynamic limit in comparison to collisionless and intermediate limits, 
where the decay is very large. At the other end in the collisionsless limit, the decay is at a 
slightly slower rate in comparison to Tc = 0.2 /m, implying that the maximal conversion rate 
is for Tc around 0.2 fm. As mentioned earlier, one peculiarity of the model is that the field 
energy is continuously created due to the recession of the plates even as the field itself decays 
to produce particles. To study this, in Fig-4.12 we present field energy per imit transverse area 
as a function of ordinary time. For Tc = 0.2 fm this energy is much less than that at = 5 




61 



Figure 4.11: Decay of the chromoelectric field as a function of proper time (in units of fermi), 
for Tc = .2/m (solid line), for = .001/m (upper dashed line) and for Tc = 5fm (lower dashed 
line). 


fm and 0.001 fm. This demonstrates that the conversion is indeed dominant at = 0.2 /m. 
This is more clearly displayed in Fig-4 13, where we have plotted the ratio of particle energy 
per unit transverse area to field energy per unit transverse area. 


In Fig-4.14 we present the evolution of temperature at Tc = 0.2 fm. The maximum temperature 
we obtain m the color flux-tube model is around 300 MeV. Finally, a brief comment on the 
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Figure 4.12: The field energy/unit transverse area(in units of GeV//m^) for = 0.2 fm (solid 
line), for Tc = 0.001 fm (upper dashed line) and for = 5 fm (lower dashed line), as a function 
of ordinary time (in fermi). 



Figure 4.13: The ratio of particle energy per unit transverse area to field energy per unit 
transverse area for = 0.2 fm (solid line), for Tc — 5 fm (upper dashed line) and for = .001 
fm (lower dashed line), as a function of ordinary time (in fermi). 
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distribution functions for quarks and gluons with explicit color dependence (Fig-4.7 and 4 8). 
As discussed earlier, it may be seen that the equilibration is fastest for 0 = 7r/2. where there is 
no background effect. At 6 = t!-/2 the thermal equilibration occurs at a common time r ~ 1 0 
fm for both quarks and gluons. 


4.4 Conclusion 

We have studied the production and the equilibration of a genuinely non-Abelian plasma with 
the color degree of freedom incorporated in both the source term and the background term 
in the transport equation. With the realistic gauge group SU(3) that we have considerd here, 
the distribution functions for quarks and gluons get defined in the extended phase space of 
dimension 14 We have added the gluonic component into the color flux-tube model, which was 
so far absent 

The role of the color degree and that of gluons was found to be substantial, and non-trivial. The 
plasma is denser than the SU(2) counter part but is comparable to the Abelian case. However, 
it IS much cooler than the abelian plasma, where the corresponding temperature is ~ SOOMeV. 
The value of the effective charge is larger than that of its counterpart m SU(2), contrary to naive 
expectations. The peak energy per particle is ~ 2 GeV which indeed is rhe demarcating scaie[37] 
between soft and hard processes. In short, the non-perturbative aspects of the evolution of QGP 
offer a rich variety of results which will have to be combined with the perturbative studies in 
order to abtain a complete description of the production and evolution of quark-gluon plasma. 



Chapter 5 


Dilepton production from a thermally 
equilibrating quark-gluon plasma 


As mentioned earler, some of the major proposed signatures to detect quark-gluon plasma 
in URHIC are 1) Jjiv supression[5], 2) electromagnetic probes such as dilepton and direct 
photon production[6, 7]. and 3) strangeness enhancement [8]. It has been suggested that for an 
equilibrated quark-gluon plasma, J jij) supression, due to screening, is a good probe[5]. However, 
due to the uncertainties in the study of production and evolution of QGP in these is experiments 
the results obtained by screening can not be taken too seriously. There are also calculations 
of J/ip suppression in the equilibrating quark-gluon plasma[9, 10] using short-distance QCD 
which are different for different models (see chapter 7). The supression of J/^. infact, has been 
observed in reactions where there is no quark-gluon plasma phase, such as in p-A collisions and 
in collisions of light nucleus [55]. This suppression is well explained by the nuclear absorption 
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of J/'0[56] Although, recently, NA50 collaboration[57] reports an excess in the supression of 
J/0, it is still not clear if a quark-gluon plasma phase has formed in this collisions. There are 
proposals that the data are explained by a deconfined partonic medium[58], or by a medium with 
high density [59]; but there are also other calculations which explain the data without assuming 
any quark-gluon plasma phase[60]. As far as the Jfxb supression in nucleus-nucleus collision is 
concerned, many aspects of it has to be studied in greater detail before unambigious conclusions 
can be drawn about the existence of quark-gluon plasma. It thus appears that there is a need 
to study more than one signature if one has to detect QGP. With that in mind we investigate 
another signature, i.e. dilepton emission. 

Dileptons and single photons have long been proposed as useful probes of the plasma[7]. Once 
produced, they hardly interact with the strongly interacting matter and thus carry the details of 
the circumstances of their production. There has been, however, uncertainties in the calculation 
in their spectra. This is due to lack of detailed knowledge about the space time evolution 
of quark-gluon plasma. Recall that the various stages by which the complete evolution of 
quark-gluon plasma is described in URHIC are, i) pre-equilibrium, ii) equilibrium, where one 
actually studies the equilibrated quark-gluon plasma, iii) cooling and iv) hadronisation. The 
pre-equilibrium stage of the collision which leads to thermal and then chemical equilibrium has 
a crucial role to play. From this point of view it is necessary to study what happens to the 
dilepton emissions from different stages of evolution of QGP, rather than estimating it in an 
equilibrated quark-gluon plasma. 

We study the dilepton production in an equilibrating QGP employing the distribution functions 
determined in the previous chapter. The dominant process which produces a dilepton pair 
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IS 


q + q i'^ + I 


(5.1) 


where 7* is the intermediate virtual photon. The dilepton emission rate in a space time volume 
d‘^x IS written as 
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Here P^ is the four momentum of the lepton pair, with Pt its transverse momentum and Mr{= 
■\flVP + Pt^) ds transverse mass. M is the invariant mass of the dilepton pair (M^ = P^P^) 
and Vrei = AP / 2 E 1 E 2 is the relative flux velocity of quark and antiquark pair in the above 
process. The distribution function of quark and antiquark is a function of the color charge Q. 
dQ is the integral of the the color part of the phcise space. The dilepton production cross section 
for the above reaction is 
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where F, = and mi is the mass of the lepton. N is the color averaging factor 

(corresponds to the volume of the compact (color) space), Ns is the spin degeneracy {Ns = 2s + l) 
and ej is the fractional charge of the flavour. For a dilepton pair of large mass the emission rate 
in the midrapidity region (Y = 0) of lepton pair is given by [61]: 
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withpt2 = - 2MTPnch{q - (i) + ph sh^2 = :^{MTshq-pnshPi) andpi = AM\\lTch{q- 

q: Pt]-K In equation (5.4) we have used = TtPpdTTdq (where R is the radius of the 
radiating region). q(= 1/137) is the coupling constant of the electromagnetic interaction. 
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It has been shovvn[62] that the thermal dilepton rate does not depend on M and Pt separately, 
but only on Mj. This is known as Mr scaling in QGP. However this Mr scaling is violated 
in an equilibiating plasma, as can be seen from equation (5.4). In this case the dilepton rate 
depends on both Mf and Pj This would also be the case if one takes the transverse expansion 
of the plasma into account [63]. 

In our calculation we have taken R = 7 fm corresponding to Au-Au collisions. We have con- 
sidered two different cases corresponding to two different relaxation times, = 5fm, and 
Tc = 0 2/rr4. The relaxation time = 5/m corresponds to collisionless limit and = 0.2/m 
corresponds to a more realistic limit of equilibration [38]. We have calculated the dilepton rate 
as a function of Mr for Pr = 0 and 1 GeV respectively. In Fig-5.1 we have presented the 
dilepton rate at = 0.2fm. As" can be seen from Fig. 5.1, the dilepton yield becomes smaller 



\tj(GeV) 

Figure 5.1: Dilepton rate from pre-equihbrium stage, as a function of Mr for =0.2 fm. 

for higher values of transverse momenta. In Fig-5.2 we have shown our results for — 5fm 
which corresponds to the collisionless limit. The rate at higher Mj is found to be larger than 
that at Tc = 0.2/m. This is because the average energy per parton is higher at = 5/m than at 
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Tc — 0.2 fm as observed in the last chapter. This is because in the collisionless limit the energy 
per parton is larger (~ 4 GeV) than that at = 0.2 fm (~ 2 GeV). 

In Fig-5 3 and Fig-5.4 we have compared our results with the corresponding Drell-Yan predic- 
tions. The Drell-Yan results are taken from the reference[28]. It can be seen from Fig. 5.3 that 



My(GeV) 

Figure 5.2: Dilepton rate from pre-equilibrium stage, as a function of Mr for =5.0 fm 
(collisionless limit). 

the dilepton rate from, pre-equilibrium stage is larger than the Drell-Yan production for very 
small transverse momentum {Pj — OGtV). This is true when Mt < 2GeV. For Mj > 2GeV, 
the Drell-Yan production dominates over the pre-equilibrium dilepton production. For a large 
transverse momentum ( Px ~ iGeV*) of dilepton pair, the production from the pre-equihbrium 
phase IS lower than the Drell-Yan production in the whole range of dilepton transverse mass. 
This can be seen from Fig-5.4. This is contrary to the earlier findings [64-66] where it is shown 
that dilepton production from the pre-equilibrium stage dominates over the Drell-Yan produc- 
tion. However the production of partons are not too hard in this model, and one expects to 
have such a low rate in the dilepton spectra than the Drell-Yan emission. 
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Figure 5.3. Drell-Yan and pre-equilibrium dilepton, rate as a function of Mr for rc=0 2 fm 
(Pt = 0 GeV). 


Summarizing the chapter, we have calculated the dilepton spectra in ultra relativistic heavy-ion 


collisions within color flux-tube model, with non-abelian features explicitly incorporated. In 


contrast to the earlier studies, we have found a smaller rate in dilepton production in the pre- 


equilibrium stage at high Mr and Pr than the Drell-Yan production. The production is larger 


in the collisionless limit of the plasma. An experimental measurement of the dilepton spectra m 


URHIC can shed light on the collision time (tc) and possibly on the initial field energy density 


that we have assumed in our model. 






Chapter 6 


A new source term for quark production 
from a space-time dependent 
chromoelectric field 


6.1 Introduction 

Over years there has been several devlopement on the production of charged particles from 
electric held, a phenomenon which was discovered nearly forty years ago by Schwinger [23]. By 
now the production of particles from a constant background electric held is extensively studied 
both theoretically and experimentally[67]. However, the subject matter of particle creation 
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from a time dependent field, and the dynamics of the ‘field-particle' system is relatively new 
and IS not fully solved. This aspect has been a subject of renewed and central interest in the 
context of heavy-ion collisions and QCD due to the suggestion that the receding nuceli might 
pioduce a strong chromoelectric field which creates quark/ anti-quark pairs and gluons[19]. As 
the particles (quark/anti-quarks and gluons) are continously produced from the background 
field, the field becomes time dependent during the evolution of the system after the two nuclei 
collide with each other. Hence the source term for the production of particles which is taken in 
chapter 3 and 4 incorporating a constant field m its calculation is not strictly applicable and a 
determination of the particle production from a time dependent electric field becomes necessary. 
Indeed a complete quantum mechanical treatment of the problem is rather too involved. For our 
purpose, it should be sufficient to derive how a space-time dependent field produce particles as 
a function of time. The dynamics of the field is of course is coupled to that of the particles, and 
may be determined in a consistent (although not complete) manner by solving the transport 
equations as we have done in the previous chapters. 

In this chapter we derive a source term for quark/anti-quark and gluon production from a time 
dependent electric field which is needed in the non-abelian transport equation[45], with the 
purpose of eventually studying the production and equilibration of quark-gluon plasma in ultra 
relativistic heavy-ion collisions. As described in the preceeding chapters the source term one 
needs is, which is the probability for particle production per unit time per unit volume 

per unit transverse momentum in the phase space. 

Before deriving a source term for a time dependent electric field there are several physical 
features that needs to be clarified. First of all Schwingers’s mechanism[23] incorporating a 
constant field in its calculation is an exact one loop non-paeturbative result, which can be 
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understood in terms of a semiclassical tunneling process[6S]. However, for a time dependent 
electric field such a calculation is not necessary and particle production can directly take place 
by the perturbative mechanism. It can be further noted that a constant field can not produce 
pairs by this perturbative mechanism. Secondly, it is not possible to obtain a source term 
usual phase space from a quantum processes This is because any attempt to 
obtain such a source term within zero uncertainty will violate quantum mechanics. Although a 
source term in the usual phase space has been proposed[68] by rewriting the Schwinger’s source 
term as an integral over pj, such proposals are ad hoc and are inconsistent with the minimum 
uncertainty requirement. In this chapter we overcome this problem by working in the coherent 
basis of coordinate, momenutm and color. 

A coherent state (ja >) is the eigen state of anihilation operator d(x,p) defined by[69] 

a|a >= ajo: > (6.1) 

with a being its eigen value the real part of which is proportional to the expectation value 
of cordmate operator and the imaginary part is proportional to the expectation value of the 
momentum operator. W.r.t this coherent state the expectation value of coordinate operator x 
and momentum operator p can be written as 

X =< Q;|f |a > (6-2) 

P =< a\p\a > (6.3) 

with AXAP = h/Arc, which satishfy the minimum uncertainity. Here {AX^ = < a\x\a >^— < 
Q|i"|a > and [AP^ = < a|p|a >^— < a|p^|a >. We have used this coherent basis (q > (= 
|,A, P >) to derive a new source term which does not violate principles of quantum mechanics. 
To deriv'e a source term for parton production from a time dependent electric field and to further 
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study its dynamics we have restricted our calculation to a finite time interval in the S matrix 
expansion (see below). 

In the next section we derive the source term for the case of a purely time dependent electric 
field. Pending the determination of E{t) from a self consistent study, we have employed a specific 
form for the electric field to get an understanding of the dynamics of the particle production. We 
have compared the results of this source term with that of Schwinger’s source term (obtained 
for a constant field) by puting the same time dependence in the field (the time dependence in the 
electric field in Schwinger’s source term was put by hand in chapter 3 and 4). In section 3 of this 
chapter we have derived a source term for quark production from a t{— dependent 

elctnc field. These source terms will be used in the transport equation to study production and 
equilibration of quark-gluon plasma in ultrarelativistic heavy-ion collisions. 


6.2 Source term for a purely time dependent electric field 

We start our derivation with a determination of the Feynmann amplitude for the production of 
quark from a time dependent electric chromoelectric field. The Feynmann amplitude is written 
either in momentum space or coordinate space. However, we need a source term (5^^^) in the 
usual phase space of coordinate and momentum from this amplitude. For this purpose, we have 
piojected the Feynmann amplitude (in momentum space) to a coherent basis of coordinate and 
momentum. .A.s described earlier this coherent state is the eigen state of anihilation operator 
and the expectation value of the coordinate and momentum with this coherent states satishfy 


the minimum uncertainty. 
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For a particle at mean position X, having mean momentum P, and color charge Q the amplitude 
in the coherent basis can be written as 

<Q,X,P\S\0>= J (fk<Q,X,P\k><k\S\0>, (6.4) 

where [.Y, P, Q > is the coherent state of cordinate, momentum and color. Here the coherent 
state IQ > for the color is defined as \Q >= ea:p(^Q“T“)|0 > (see Eqn 6.14). To the leading 
order in perturbation series we have 

5 = ^ dt'd^xjCx{t',x). (6.5) 

with /lj(i T x) being the interaction lagrangian density. The time [t = 0) refers to the initial 
state whmc there is no production (this replaces t = — oo in the usual scattering processes). 
For our dynamics here we have limitted the time integral to a finite limit t and not to co. For 
quark/anti-quark pair production the interaction lagrangian density is 

£r = (6 6) 

with i/q b(ung the Dirac field for quark and being the external color field. In the above 
equation T" are the generators with a = l..,.,8 and = 1,2,3. For quark/anti-quark pair 
production the corresponding amplitude is 

< Qi, -Vi , A'2, Pi - PilPjO >= Sjj J d^kid?k 2 < Qi, Q2, ^1, AI2, Pi, P 2 |fci, ^ 2 , J > 

< ■31,52, All, A:2|S'|0 >, (6.7) 


where ,Yi, AQ are coordinates Pi, P 2 momenta and Si and S 2 are the spin indices of the quark 
and anti-quark respectively. Here < i,j,SuS2,ki,k2\S\0 > is the usual Feynmann amplitude for 
quark/anti-quark pair production with i and j being the color index of quark and anti-quark. 
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After summing over spin and color the absolute square of this amplitude is written as: 


|< QuQ2,XuX,,P,,P,\S\Q > 1 ' = J d%d^k[d^hd% < QuQ2,XuX2,PuP2\kuk2,i/X, 

< AA, Al2,Pi,P2 >*< OlPIfc^, A;2,Si,S2,/,m >*< z,j/,Si,S2,^i,^2|*S'|0 > 


In our calculation, we ha\e made a suitable gauge choice such that = 0. With this 

choice the potential in the the above equation is related to the electric field by the simple 

expression. 


E^{t) = - 


dt 


(6.9) 


Substituting 

.4(t,x) = j d^KA{t^K)ex-p(—iK-x) 


( 6 . 10 ) 


we obtain (for massless quark) 


^51,32 1^1) ^2 ^ < ^1 1 ^2|5'|0 > — /(ti, fci, A:2, A:2)) (d.ll) 

where /(ij , Aii , A; 2 , /;(, A;!) is given in the appendix for a fixed value of a in . The expression 
foi roller sum (C(Q)) is not included in /(ti, ^ 2 ) Aii, A: 2 , A:(, A:^) and is presented separately (see 
Eqn (6.14)). With the above founer transformation (Eqn (6.10)) the time dependence in A 
appears finally in the source term. In most of the earlier studies[71] the source term is derived 
witii some special form of the electric field. Such source terms are not necessarily applicable 
m the transport equations to study the production and equilibtation of quark-gluon plasma. In 
our study, here no specific form for electric field is assumed. 

For the case of a single quark we project Feynmann amplitude to a coherent basis \X, P,Q > 
rather than projecting it to 1 Ai, A: 2 , Pi, P 2 , QuQ 2 >■ The later is needed if we are interested in 
pair (quark/anti-quark) production and corelations. For quark production at mean position X 
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having mean momentum P and color charge Q we have 


1<Q,.Y,P|5|0>|' = 


k <Q,X,P,\k,,.' ixj,. k[\Q,X,P>‘ 


< 0\S\k[,kJ,l>^ < l,i,k,,k\s \0 > . ( 6 . 12 ) 


The above equation can be written as 


l<Q,X7^|.hlO>l' 


k<Q\t><j\Q < X, P, |^-l > < k[\X,P > 


< /riX|5“|0 > < >('6.13) 


As was done in the previous chapters we fix the direction of the electric field and we take it to 


b(' along 8 dinvtion in the color space (a =8 and b= 8 ). The color factor C(Q) becomes, 

C{Q) = < (5b >< j\Q = ^[1 + 3| < gi3 > |'], (6.14) 


with \Q >= (ixp{iQ°-T'^)\Q > 


Following the reference[70l we arrive at the expression < (5|3 > = cosd ex-p{i'y)-, so that C{Q) = 
;^[1 + 3 ras'^/?]. This is the aniverasl form of the color dependence independent of the nature of 
the fleet ric Jicld. Here 0 < 9 < it 12 and 0 < 7 < tt. 

For a purely Lime dependent potential A{t.ki + k 2 = 0) we obtain 

1 < Q, X P|S|0 >1^ = CiQ) J d^k,d^k[2d^k < X P; \ki > < k[\X, P >’ 

[ dti f dt2f{tut2,ki,k,k[,k)exp{i{k°-Tk°)ti)exp{-i{k[° + k°)t2)S^{ki+k)5^{k[+kJ.6.15) 

Ji} Jo 

I'he colKM-ent state wave function < X, P\k > which satisfies minimum uncertainty is 


< X, P\k >= N exp{-{P - kf l2Q'^)exp{ik ■ X), 


(6.16) 
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with N being the normalisation constant. Putting the the above form for the coherent state m 
eciuation (6.15) we get 


l< 0. A'.PISIO >1= = ,V=C(Q) I <ekexp{-ik-P)yQ^) j' dt, £ di, 

- t,))~[k\A{t,) ■ A(h)) - {k A(t2)){k ■ A{h))], (6 17) 

With k — |A,| (for massless quark) Without any loss of generality one can assume z-axis to be 
along the direction of A and obtains 


|< ()), A, P|b]0 >1 — C{Q)N^ j (Pkex'p(—{k — Py/Q'^) dti J dt2{A{ti)A{t2))ex'p{2i\k\{ti — t2) 


-[(1 - cos^e^ 


/f: 


wit.h Of, being tlic angle between z-axis and k. This is the probability of quark production per 


unit volume per unit momentum having a color charge Q from a time dependent electric field 
E{t) related to A{t) by Eqn (S^j^The required source term is then 


dN 


dtdA 


iiv _ r 

X(PBr ~ J 


dP, 


d\ < (5,X,P|5|0> 


dt 


(6.19) 


I'u make sun' t.hat a constant potential and a constant field (which is obtained from a linear time 
dependent potential (Eqn. (6.9)) do not contribute to the above time dependent source term 
w(> project out only the time dependent part by the replacement A{t) A{t) — A(0) — tA'{0) 
in the Eqn. (6.19). .-Vfter singling out the contributions from this constant potential and the 
constant held we obt ain a genuinely time dependent source term from Eqn. (6.19). 

Now by putting some specific form of the electric field we compare the results obtained by this 
sourc(“ t.erm with that of Schwinger’s source term. For the choice A(t) = exp(— t/to) we present 
the results of both the cases below. W^e have choosen Qx and Qs iti the coherent state wave 
function (Ecin. 6.16) corresponding to a lAl dimensional expanding plasma. For Qx the scale 
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IS the transverse diameter of the nucleus 14 fm for Au nucleus) and for it is determined 
from the total length traxeled by the two nuclei (2z or 2t, as nuclei travel almost at a speed of 
light) In Fig-6.1 we have plotted the production rate for 0.5 GeV as a function 

of i The production rate increases as time increases in contrast to the Schwinger’s results 



P’lgiire fi.l; The rate of quark production obtained by the new source term 

from <i tiiiK' dependent electric field for Pt = 0.5 GeV, to = 0.5 fm and 9 = 0 

whuh IS shown in Fig-6.2 for the same input parameters. This is because, in the newsource 
term we have taken a finite limit in the time integration in the feynmann amplitude. As time 
increases, t.he ab.solute value of this amplitude increases and hence a large production rate is 
ol)t<uaed. However, the schwinger formula is valid for a constant field. But for a comparison 
with our new source term we have taken the same time dependence in the electric field in both 
the cas('s. In Fig-6.3 we present the results for Pt = 1.0 GeV. The production rate decreases 
as Pt increases, which is in accordance with the Schwinger’s result (presented in Fig-6.4). 
I'liis IS more clearely seen in Fig-6.5 where where we have plotted the production rate of the 
new source term for Pt ” 5.0 GeV. To get a feeling on the total production of particles rather 
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*) 2: rate of quark production obtained from Schwinger’s source 

U'lru vvit.l) the same time dependent electric field (as described in the text) for Pt = 0.5 GeV, 
tfj = 0.5 fm and 0 = 0 

than it.s i)ro<lu(ti(;n rate we have plotted the total production In Fig-6.6 and 6.7 for 

/’;• := 0.5 a,nd 5.0 (le.V up to t = 2.5 fm. It is seen that total number of particle production per 
unit voluiiH' per unit transverse momentum is an increasing function of time. This has a strong 
miiiact on ih<‘ pat tide production in relativistic heavy-ion collisions where the plasma would be 
very dtuise eviui at a later time (up to 2.0 fm). Such a calculation of the quark-gluon plasma 
pioduition and it.s evolution within the frame work of non-Abelian transport equation will be 
[lerformed and reported elsewhere. 
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iMgurc 1 lie rate of quark production (GeV"^)) obtained by the new source term 

Frorii a time dependent electric field for Pt = 1.0 GeV, to = 0.5 fm and 6* = 0 



Figure 6.4: The rate of quark production obtained from Schwinger’s source 

term with the same time dependent electric field (as described in the text) for Pt = 1.0 GeV, 
to = 0.5 fin and 8 = 0 
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{•'iguic f) .") riic rate of (juark production ( ild^^Xd^Pr obtained by the new source term 

fioin .1 tiriH* d(‘[)(‘udoiit (doctric field for = 5 0 GeV, to = 0.5 fm and ^ = 0 



Fif>;ur(‘ ti.B: I'lu' tjuark production )) obtained by the new source term from a 

time d<‘[>riu!<!nt tdoctric field for Pp = 0.5 GeV, to = 0.5 fm and 5 — 0 
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Figure ().7: The quark production {-jzj^iGeV )) obtained by the new source term from a 
time dependent electric field for Pt = 5.0 GeV, to = 0.5 fm and ^ = 0 
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6.3 Source term for a r dependent electric field 


Ah seen in the previous chapters, for a 1+1 dimensional expanding quark-gluon plasma we 
have employed Bjorken’s boost invariance formalism according to which the longitudinal boost 
arc the symmetry operations on the particle production. The boost invariant parameters are 
r = Pt and ^ = r] ~ y with t] = tanh~^{z jt) and y — tanh~^{P3l Po)- The source 

term for particle production (^, Px)), needed m the transport equation depends only 

on these boost invariant parameters. In this section we derive such a source term which is 
needed in the transport equation to study production and equilibration of quark-gluon plasma 
in ultrarclativistic heavy-ion collisions. 


In terms of these variables we get dtdz =■ r dr drj and 

S = [ T'dT'dr]dxdyCx{r\r],x,y). 
Jo 


(6 20 ) 


s r 


(this equation is analogous to Eqn. {P) for a purely time dependent electric field). To be 
consistent with the previous chapters we work in the Lorentz gauge d^A>^ = 0 with = 




duG{r) 


= 0,3. Due to this choice the chromo-electric field depends only on the 


boost invariant quantity t which is given by 


E{r) = [^,+li]G{r) = A + hff(r). 


Following the same procedure as was 


dr^ ' rdr^''^ ' ^dr 

done for the time dependent case above, we get 


( 6 . 21 ) 


l<g,x,P]5lo> 


= s,,.; J d%d%d^k <Q\^><J\Q 1^1 > < ^ > 

r dnn r dT2ndmdr]29{ku k, k[, k,r,,T2, Vum) exp(Kfc? + k^)TiCoshrj,) 
Jo Jo 

exp(-i(A:(° + k°)r2Coshr,2)exv{-i{ku + k,)Tisinhyi)exp{i{k[, + K)T2sinhri2) 

5\k'^,T + kT)5^{k[T + 
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wheie the expression for g{ki,k, k[,k^ Ti, r2, 771, 772) is given in the appendix for massless quark. In 
the above equation we have used the fourier transformation A{t, z, Xt) = f cPKtA{t, z, Kt) exp(-j AV 
Xt)i with Aj = 0) for the vector potential which depends only on t and z. Eqn (6.22) 

can be further simplified and is written as 

1< Q,X, AIAIO >1' = N^C{Q) J dk,,dk[J%dhexp{-ikT - Pt)VQI) 

exp(-((A;i2 - + [ku - PzY)l‘2.Ql) driTi dr2r2dr]idr]2 

g{ki,k,k[,k,Ti,T2,r)i,T]2)exY>{i{k° + k°)TiCoshrii) ex^{-i{k[° + k°)r2Coshg2) 

exp(— 7(^12 + k:;)Tisznhgi)ex^{t{k[^ + kz)T2Sinhg2) ex-^{—i{ku — ^u)Al3), (6.23) 

with IP = \l[k‘^ + kl), k\ = \J{k‘^ + klP) and kf = \J{k‘^ + k^P). Substituting Xz — rsinhg and 
Ai = Pr^'^inhy with = 77 — y we get 

|< Q,X,P\S\0 >f = N^C{Q) j dkudk[J^kTdk,exp{-{kT - Pt)VQt) 

exp{-{{ku- PTSinhyY + {ku- PTSinhyf)/ 2 Ql) [ dTin f dr2T2dgidrj2 

Jo JO 

Sikuk,k[,k,gi,g 2 )G'{Ti)G'{T 2 )ex^{i{k° + k'^)TiCoshgi)exp{-z{k'° + k^)r 2 Coshrj 2 ) 
exp(-i(A'i;: + k^)rxsinhgi)ex-^{i{k[^ + K)T2S%nhg2)&yA>{-J^iz - k[J}Tsinh{^ + y)) (6 24) 

where 

S{ki, k, A:, 771, {ti)G' {T2) = y(A:i, fc, Aii, A;, ri, r2, 771, 772). ( 6 . 25 ) 

Integrating this expression over P3 and taking derivative w.r.t. r we arrive at the new source 
trem, 

>l^ (6.26) 

dtdP'xdPPT dr J 

which depends only on the boost invariant quantities r, ^ and Pt- This source term will be 
used in the transport equation to study the production and equilibration of quark-gluon plasma 
m ultrarelativistic heavy-ion collisions. This work will be taken up separately. 
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6.4 APPENDIX 


For m = 0 one obtains 


f{ti, t2,ki,k2,k[,k!2) — + I2) 

JD 


(6 27) 


where 


B = ^/k°k[%^k!2° 


(6.28) 


and 


= [{k, ■ A){k2 ■ A-)i4k[^k'^° - {k[ k2)) + (fci • A)(k2 ■ A*){k, ■ k'^) + (^1 A){k[ • A*)(/j( • ^ 2 ) 
-{k\-A){k'^ A-){kr-k2)]A[{h-A*){k2 A){Ak[%'^^ -{k[-k2))^-{k[-A*){k2-A){k, k'^) + 
{ki ■ A'‘){k\ ■ A){k[ ■ /C 2 ) — {k[ ■ A*){k2 ■ A){ki ■ k2)] + {A ■ A’)[4A;( k2 {k\ ■ ^ 2 ) + (^1 ■ ^i)(^2 • ^ 2 ) 


+(*, . fe)(4; . i-i] - (k, . ki)(Li . 42)1 + (((*1 • ■ ^') - (^1 ■ ■ K) + ({k2 


( 4 ; ■ A) - (fcj ■ A)(A:; . A-))(4i 4^.29) 


and 


I, = (4, ■ A)(42 . A-)[24;"4'° - ( 4 ; . 4ai + (4. A)(4' ■ A-)[(42 • k[) - 2k[X°] + (4i ■ A')(42 • A) 
[24:°4'° - ( 4 ; . 4 ;)] + (4, . A-)(4; • A)[(42 ■ 41) - 24;°42“1 + (4; • A){4i ■ A-)(42 ■ k',) - 

(^k[ ' A*){ki • A){k 2 • ^ 2 ) (^1 * ^ * ^ 2 ) ^2 ] ^ )(^2 * * ^ 2 ) 
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+2k,%2°] + {k[ • A)(/c2 ■ A*)[{k, ■ k',) - + {k[ • A)(fc' • A*)[-(A:i • Ars) + 2ky% 


+{k2 ■ A)ik'^ ■ A‘)pi • k[) - 2k,%[°] + (fc; • A){k2 ■ A-)[-(fci k[) + -A-A* 

[2k,°k'^\k2 ■ k[) - 2k,°k2°{k[ ■ k'^) + fei°fc2°(A:x • k' 2 ) - k[°k'^°{k, • k{^.30) 


More generally /(ii, i 2 , ^ 1 , A: 2 , /cj, A: 2 , A:^) where A = A°’andA* = A*”. 


and 


g{ki,k[,k2,n,T2,r]i,rj2) = — 


( 6 . 31 ) 


with 

C = k°^/k°k[° 


( 6 . 32 ) 


and 

/ = A ■ A*[2A:f (All ■ Ai() + 2Ax°A2°(it2 • A:() - Qk[°k2°{h k2)]+{k, ■ A){k2 • A*)(6A:(°A:2°) 

+(A;2 • A)(A;i • A"‘)(6A:(°A:2°)— (A:( • A)(A;i ■ A*)(2A:2°A:2°) + {k[ ■ A‘){ki ■ A)(2A:2°A;2°) 

-{k[ ■ A)(A:2 • A*)(2A:i%°) - (A;( A’)(A:2 • A){2ki°k2°) + A°[(A:i • A*)2k°{k2 ■ k[) 

+{k2 A*)2kl{h k[)-{k[-A*)2kl{h-k2)]AA°’^[{k[-A)2kl{k2-ki) 

-{ki • A)2k°ik2 ■ k[) - {k2 ■ A)2kl{k2 ■ A:()] (6.33) 

( if tXf ■ i I ' ] 

Here Aq = -sinhr,^G'{Ti), A 3 = coshg^G' { ti) , Aq = -sinhr] 2 G'{r 2 ), A^ = coshr) 2 G' {T 2 ) and 
K ‘ 

other A’s are zero. 



Chapter 7 


suppressions in a thermally 
equilibrating quark-gluon plasma at 

RHIC 


7.1 INTRODUCTION 

As mentioned earlier, dileptons and single photons have long been proposed as useful probes 
of the plasnia[7], as once produced, they hardly interact with the strongly interacting matter 
and thus carry the details of the circumstances of their production. However, as they are also 
produced via hadronic decays in the later stages of evolution, where plasma expands and cools. 
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it becomes difficult to distinguish between them. In that sense the production of cc is a clean 
process, as it is not produced m the later stages due to its heavy mass. The production of cc 
takes place mainly at the hard vertex. Once produced, cc will evolve in to open D mesons and 
charmonium states, such as J/4>, x^ while travelling through this dense phase of matter. 
Unlike hadronic collisions, where an enhancement of J/0 is found at tevatron energy[72], both 
cc and J jil) get suppressed while traveling through different stages of quark-gluon plasma. 
There might be thermal charm productions in the very early stage of the plasma, where the 
temperature is very high, but no production of cc can occur in the later stages due to decrease 
m the temperature. 

We again mention that, the various stages by which the complete evolution of quark-gluon 
plasma is described in URHIC are, i) pre-equilibrium, ii) equilibrium, where one actually studies 
the thermalised quark-gluon plasma, iii) cooling and iv) hadronisation. The production of an 
equilibrated quark-gluon plasma in stage (ii) crucially depends on the pre-equilibrium evolution, 
i.e. on stage (i). Matsui and Satz[5] have suggested the suppression of in the equilibrium 
phase. In their study, the debye screening length which is calculated from lattice QCD, is found 
to be less than the J/(/> radius. This forbids the binding of cc to J/i/^ in the equilibrated quark- 
gluon plasma. However, such a study is not available for an equilibrating quark-gluon plasma. 
There could be many interesting effects in the pre-equilibrium stage, when cc or J/i/^ travel 
before reaching the equilibrium stage. As high energy deconfined partons are also present 
in the pre-equilibrium stage, J/^ is suppressed due to its interaction with these deconfined 
partons. Since the J/U dissociation cross section inside a deconfined partons is very different 
from that inside a hadronic gas (see section II), the study of J/V> suppression can then provide 
us information on color deconfinement in the quark-gluon plasma. 



91 


Another outcome of the pre-equilibrium study is the prediction of the equilibration time, the time 
at which quark-gluon plasma equilibrates. A detailed knowledge of equilibration time and J jil) 
formation time is very crucial in determining the Jjip suppressions in URHIC. This is because, 
a large equilibration time may permit & J jip formation from cc pairs before plasma equilibrates. 
Then the interaction of a fully formed J/0 with deconfined partons is prefered than screening. 

If the J/'i/) formation time is greater than the equilibration time or comparable to it, then cc 
interacts with the deconfined partons before screening starts operating in the equilibrium stage. 
The interaction of a cc with the partons before plasma equilibriates is not studied yet, and such 
a calculation is beyond the scope of this thesis. 

.As mentioned earlier, a careful study of pre-equilibrium evolution of quark-gluon plasma is the 
consistent way to study production and equilibration of quark-gluon plasma in URHIC. From 
this point of view it is necessary to study what happens actually to the J/i/j suppressions at 
different stages of evolution of QGP, rather than estimating it in an equilibrated quark-gluon 
plasma. We study here the suppression of a fully formed J/V’ in equilibrating quark-gluon 
plasma using short distance QCD. The J/i/> may also be suppressed in later stages due to 
inelastic collisions with hadrons, but we do not consider it here. This suppressions will be less 
due to lesser value of the hadron- J/i/) inelastic cross section (described in section II). 

However, J/^ suppression is also measured m p-A collisions where there is no existence of 
quark-gluon plasma[73, 55]. At these experiments the suppression of J/i/^ is due to the presence 
of a nuclear medium. The prominent mechanism by which these datas are explained at SPS is 
via nuclear absorption, the suppression of a nascent J/i/; before it forms a physical resonance[56, 
74]. Within this mechanism, one does not assume the existence of a quark-gluon plasma phase 
or a highly densed deconfined partonic medium. This is the case with some light nuclei collisions 
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at SPS[73, 55]. However the recent measurements of Jjij) suppression by NA50 collaboration[57] 
(Pb-Pb collisions at = llGeV) yeild an excess suppression of which is not explained by 
the above conventional approach. Presence of a nuclear medium may not be enough to explain 
this data. There are speculations that this suppression is due to the e.xistence of a deconfined 
partonic medium[58] or a high energy density matter[59]. 

Whether an equlibrated plasma has formed in Pb-Pb collisons at SPS(-y/i = llGeV) or not, is 
debatable, but at RHIC(Au-Au collisions at = 200Gey ) we may go closer to the equilibrated 
quark-gluon plasma. In this chapter we estimate the J/V’ suppressions at RHIC for a thermally 
equilibrating quark-gluon plasma. 

However before studying charmonium suppressions in quark-gluon plasma it is essential to study 
the corresponding situations in hadronic collisions. The experimental analysis of quarkonia pro- 
duction in high energy hadronic collisions at tevatron energy[72] revealed a drastic disagreement 
with the predictions from color singlet model[75]. A substantial fraction of quarkonium produc- 
tion at very high energy arises from intermediate color octet fluctuation during the evolution 
into final state color singlet QQ pair. Now a rigorous theoretical foundation has been achived 
by Bodwin, Braaten, and Leepage[76], who devloped a QCD formalism by marrying pQCD 
and an effective field theory within non-relativistic QCD(NRQCD). This formalism accounts 
for the production of both color singlet and color octet cc states, that evolve in to final color 
singlet quarkonium. Within NRQCD, the production amplitude is expanded in powers of both 
the strong coupling a, and the velocity of heavy quark(u" = .23 for cc system), which includes 
higher fock-state components. The long distance non-pertubative matrix element is either fitted 
from experiments or taken, in principle, from lattice QCD[77]. This technique has been used 
to study many heavy quarkonium production processes in hadronic collisons. However a space 
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time evolution of quarkonium using this technique has not been achived. This is crucial for 
stud} ing the evolution of in nuclear collisions. The existence of background chromoelectric 
rield(45] and high dense partonic medium complicates the study of the evolution of cc to J/^, 
both at RHIC and LHC. As cc takes 1.0 fm (time in cc rest frame) to form a Jfip bound state, 
it interacts with the nearby light partons and with the background chromoelectric field during 
this period. No attempts is made on this line to study production and its subsequent evol- 
ution in nucleus-nucleus collisions. This might be the consistent way to study J/V’ suppression 
in URHIC. 

Since we study the suppression of Jjtj) in equilibrating quark-gluon plasma, it is relevant to 
discuss the process of equilibration of QGP in URHIC, which is not completely studied yet. The 
rate of equilibration is different for diiferent models (see chapter 8). Some relevant models that 
describe the equilibration of QGP are parton cascade model(PCM)[17], heavy-ion jet interaction 
generator(HIJING)[78] and color flux-tube model. PCM and HIJING are both perturbative 
QCD based models, and color flux-tube model is a field theoritical model. The former models 
describe the evolution of hard and semihard partons and the later describes the production of 
soft partons vua non-perturbative Schwinger mechanism. For a complete study of production 
and equilibration of quark-gluon plasma, color flux-tube model may be combined with the 
pQCD based models. We use here the distribution of partons from parton cascade model to 
study J/'0 suppressions in a thermally equilibrating quark-gluon plasma at RHIC. A study of 
J/i/) suppressions in a chemically equilibrating quark-gluon plasma within HIJING model can be 
found in the ref. [10]. As mentioned in chapter 4 the J/i/i suppressions in a thermally equilibrating 
quark-gluon plasma within color flux-tube model will be calculated once hard partons are added 
to it via a minijet source term. Within parton cascade model, a relativistic transport equation 
is solved to study the distribution of partons by considering the direct collision processes 2 -)■ 2, 
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the inelastic collision processes 2-^1, and the decay processes 1 -> 2[17]. More explicitely, 
this is written as 

= C',(x,p) C^-^) 

where f^{x,p) is the distribution of a particular type of parton i m the usual phase space. The 
collisions term C,(a;,p) is derived from pQCD taking all the processes mentioned above. The 
initial condition on ft{x,p) is the distribution of partons inside the nucleus before the two nuclei 
collide with each other [38]. However for a complete study, the addition of soft partons and the 
effect of background field has to be taken into account. These backgruond electric field may 
not exit in the very early stage where pQCD is applicable[41], but in later times this will have 
an important role in the equilibration of the plasma[41, 40]. Such calculations may be available 
once the color flux-tube model is combined to the pQCD based models. 

In section 11 we describe the short-distance QCD and the dissociation oi J/^. Results and 
discussions are presented in section III, for RHIC energy. We summarise and conclude the 

main results in section IV 


7.2 J/i/' dissociation by deconfined gluons 

In the framework developed by Bhanot aad Peskinpg, 80], interaotioos between light hadrons 
and deeply bound, heavy quarkoninm states, such as the W, axe mediated by short-range 
color dipole interactions. For sufficiently heavy quarks, the dissociation of quarkoninm states 
by interaction with Ught hadrons is fully accounted for by short-distance QCDpg, 80). Because 
of its small size, a heavy quarkoninm can probe the short distance properties of light hadron. A 
parton based calculation of the J/t/i-hadron cross section is thus possible via an operator product 
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expansion method, simillar to that used in deeply inelastic lepton- hadron scatterings [79-82]. 
These pertubative calculations become valid when the space and time scale associated with 
the quarkonium state, rq and tg, are small in comparison to the nonperturbative scale AQg£), 
which is the characteristic size of the light hadrons, i.e. rq «< ^qcD’ ^qcd- 

For J ground state rjj^ ~ 0.2/m = {lGeV)~^ and Ej/^,{2Md — M.^i) ~ 0.64 GeV. With 
Aqcd = Q.2GeV , the above inequalities seem to be well satisfied[54] and one expect that the 
dissociation of J/?/> by hadron will be governed by the J/i/’-hadron break-up cross section as 
calculated in short distance QCD. The operator product expansion allows one to express the 


hadron- J/i/ inelastic cross section in terms of the convolution of the gluon- J/i/ dissociation 
cross section with the gluon distribution inside the hadron. The gluon- J/^/i dissociation cross 
section is given by[79, 80, 82]; 


27r 


, IStt, 


) - y(32/3) ( 2^2 (g0/e^)5 ■ 


(7.2) 


Here is the coupling between gluon and charm quark, mg the charm quark mass, and q° the 
gluon energy in the J/i/> rest frame. The heavy quarkonium is coulomb like, and one writes. 


eo = 82]. 


The gluons which are soft inside a pion are not capable of dissociating a charmonium. Hence 
they give a low tt-J/i/ cross-section. However, deconfined gluons, which carry enough energies, 
are sufficient to break a charmonium [54]. These high energy partons are present in LRHIC, 
as observed in several calculations, such as[45, 17]. These deconfined gluons can break a fully 
formed J li> that exit inside QGP. This conclusion does not seem to be affected substantially 

by uoxiperturbative effects [83]. 

We use the above g-J/tA dissociation cross section to study the suppressions in a thermaily 
equilibrating quark-gluon plasma at RHIC. For the central collisions, the survival probab- 
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ility is: 


Sg^j/^iPr) = exp[- J drugir) < Vreicr{k • v) >fc]. (7.3) 

Heie ri^(r) is the gluon number density which evolves as system expands longitudinally. Accord-* 
ing to Bjorken scenerio[2], this number density is a function of the boost invariant parameter r 
^ V ^‘^)) The thermal average g-Jl^p cross section, < Vreio'{k • u) >/., is 




(7.4) 


Here v{'= {Mx, PT-,0)j Mj/^) is the four-velocity of J/tp in central rapidity region and k is the 
four-momentum of gluon in the parton gas. In the thermal average cross-section we use the dis- 
tribution function, /(fc°,T(r)) = exp(fc°7 ' 7(T))-i gluon. Here a(r) captures the deviation from 
equilibrium, which m principle, is determined from the relation, ng{T) = f dr(p^u^)f(p,T(T)), 
by knowing the evolution of temperature and Ug , for an equilibrating quark-gluon plasma. Here 
dr = - " 3 ^- with 7 = 2 X 8(the product of spin and color degeneracy), and u( = (t/r, 0, 0, 
z/r)) is the flow velocity. However, as a(r) cancels both from numerator and denominator 
of equation(7.4), we need not know its form. What matters is ng{T) and T(r) whose evolu- 
tions are taken from PCM[17]. Within PCM, ng[T) ~ 0.72n(r), with n(r) = 


and T(t) = %0MeV{ 


X '1-0.30 

To/ ’ 


for a thermally equilibrating quark-gluon plasma at RHIC. Here 


To = .05/m. Using this distributions the integration in equation (7.3) is performed numerically 
to study the //■^ survival probability. 


7.3 Results and discussions 

In Fig-7.1 we present the thermal averaged gluon-J/^ cross section as a function of J/tp trans- 
verse momentum Pt-, for different values of temperature. It can be seen that this value is de- 
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creased as the temperature of the plasma is increased. Also for J/V’ with high Pt this thermal 



0 1 2 3 4 5 6 7 8 9 10 


P. (GeV) 

Figure 7.1: The thermal averaged gluon- J/^ cross section < Vrei^' > {mh) as a function of 
transverse momentum Pr at different temperatures. Solid line refers to T= 0.2 GeV, upper 
broken line corresponds to T= 0.4 GeV and lower broken line refers to T=0.8 GeV. 

cross section becomes less. This has a direct impact on Jj ip suppression. As Pj becomes large 
the survival probability must increase. However, due to large initial parton density formed in 
the initial stage of the plasma, these smaller thermal cross sections do not have a large effects 
on the total survival probability. This is shown in Fig-7.2 where we have plotted the survival 
probability as a function of Pt for a thermally equilibrating QGP at RHIC within PCM. 

As can be seen, the suppression is almost 100 percent for J/ip with low transverse momentum. 
Even for J/ip with high Pr(— lOGeV) the suppression is about 95 percent. This huge suppres- 
sion IS due to the presence of a highly densed deconfined partonic medium at RHIC. As can 
be seen from equation(7.3) and Fig-7.1 this suppression will be more for a plasma with larger 
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Figure 7.2: The survival probability of J/i/> in a thermally equilibrating plasma at RHIC. 


number density and lower temperature. A denser and cooler plasma is obtained once the soft 


partons studied within the color flux-tube model in earlier chapters are added to pQCD based 


models such as PCM. The meaisurement of such huge suppression will reveal the existence of 


such a deconfined partonic medium, possibly a thermalised quark-gluon plasma, at RHIC. 


However, what one actually measures experimentally is the total survival probability This 


implies that other sources of J/ijj suppression, without the existence of a QGP phase, have to 


be singled out from the total survival probability in order to see the suppressions only from this 


deconfined partonic medium. If the estimated suppression due to all sources without QGP phase 


is much less than the experimental findings, one may then hope that the rest of the suppressions 


is due to the existence of a highly densed deconfined partons. The experimental measurement 


of J/i/) at RHIC will be able to reveal all the possibilities and may explain the existence of such 


a deconfined partonic medium, possibly a thermalised quark-gluon plasma. 


On the other hand enhancement can occur due to an increase in the (thermal) charms 
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prociucod at a very high temperature, especially in the initial stage of the evolution of the 
system[bo, b6] We recall that screening of Jjxj} was originally[5] proposed without taking 
firoducLion of thermal charms into account, which was argued to be negligible at T ~ 200 - 300 
McV However, the system may approach a temperature which is as large as 900 MeV initially, 
and can therefore produce a large number of thermal charms at RHIC So to make resonable 
estimates at RHIC and LHC, one has to take into account both thermal and hard J/ip. For 
simplicity, we have not considered here those J/i/i which are created from thermal charms. 


7.4 Conclusion 

Applying short distance QCD mechanism and taking parton cascade model for the evolution of 
quark-gluon plasma at RHIC we have found a huge suppression (almost 100 percent) of J/r/i 
which are produced from primary charm quarks. In principle, the total survival probability is 
written as 

St = Sg^J/^Sother (1-5) 

where first one {Sg-j/^) is due to g-J/i/’ dissociation, and the last one{Soiher) is due to any other 
possible sources of suppression[ll-14] such as nuclear absorption, final state comover scatterings 
etc. As the hadron- J/'0 cross section is smaller than g-J/i/i cross section[54], we expect a lesser 
suppression of J/'^ in the hadromsation stage. The prominent suppression seems to be due to 
the nuclear absorption (60 percent at RHIC[9]) and J 1'^ dissociation in the deconfined. partonic 
system. Hence the measurement of a huge J/i/i suppression at RHIC may suggest the existence 
of a partonic plasma. However, many refinements, such as space time evolution of quarkonium 
production in nucleus-nucleus collision using NRQCD has to be worked out, before unambigous 
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conclusions can be drawn for the pre-resonance suppression oi J jtj). Also, the contribution from 
thermal which may not be neglible at RHIC and LHC, haa to be taken into account. It is 
only after incorporating the above features the certainty of screening can be justified. 

In any case, we have shown here the importance of the thermally equlibrating quark-gluon 
plasma on Jlijj suppression. This suppression depends on variuos models which describes the 
thermal evolution. From the mea.surement of this suppression, the validity of these models for 
thermal evolution of quark-gluon plasma at RHIC can be checked. 



Chapter 8 


A brief description of other Models and 
approaches 


8.1 Dynamical Models in heavy-ion collisions 

8.1.1 Hevy-ion jet interaction generator (HIJING) 

Heavy-ion jet interaction generator model[37] is a pQCD inspired model where the non-perturbative 
part is modeled phenomenologically The evolution of hard and semihard partons are described 
by pQCD. The intitial condition is estimated by HIJING monte carlo code in which initial dir- 
ect parton scatterings are taken into account. The produced parton gas in the central rapidity 
region approaches thermal and chemical equilibrium through elastic scatterings and induced ra- 
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diations. In an overlap region of Ifm thickness (of the two nuclei) a local isotropy in momentum 
distributions assumed and is estimated to be achieved in a time period of 0.5-0. 7 fm through 
fiv'f'-streaming evolution. It is further assumed that this is the actual kinetic equilibration ( or 
th(u-malisation) time for the partonic system. Though this assumption is not well founded, this 
is the best tliat the model incorporates in the absence of a complete calculation of the relaxa- 
tion time. Within this model the scattering processes which maintain thermal equilibrium are 
gg ^ gg^ gq gq etc.. Further chemical equilibration of partons is studied by by solving master 
rate eipiations, where gg ^ ggg and gg ^ qq processes are taken into account. Asssuming that 
the parton scatterings are sufficiently rapid so as to mainatin thermal equilibrium, and therefore 
neglecting the effects of viscosity due to elastic and inelastic scatterings, one determines the 
evolution of the energy density from the hydrodynamic equation. The master equations are 
solved self consistently[37] to obtain the time evolution of the temperature and fugacities. It is 
found that the parton gas cools considerably faster than the Bjorkens scaling law (TV = cons.), 


because the production of additional partons during the approach to chemical equilibrium con- 
sumes an appreciable amount of energy. It is shown within this model that the equilibration 
can hardly reach before the temperature drops below 200MeV in the short period time of 1-2 
fm at RHIC energy. However, at LHC energy the parton gas, especially the gluon component, 
becomes very close to equilibrium and the plasma may exist in a deconfined phase for as long 
as 4-5 fm. One observation of this model is that the quarks never reach chemical equilibrium at 


both energies. 
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8-1.2 Parton Cascade Model (PCM) 

1 lid parton cascade model[17] is a relativistic kitetic model for high energy nuclear collisions 
which is inspiied by the parton picture of hadronic collisions. The dynamics of the space time 
(wolution of quarks and gluons is studied by the interaction of partons in the framework of 
p<,‘rturbativc QCD. The time devlopment of parton cascades is taken in six-dimensional phase 
space of coordinate and momenta according to a semiclassical, relativistic transport equation 
witii a collision kernel that contains the asymptotically dominant perturbative QCD interaction 
processes among the partons. The time devlopement of the parton distributions, the number 
and type of interactions is determined by the dynamics itself. The transport equation which is 
solved is 

p‘S,F.{r,p) = SJ<«:)(r,p) (8.1) 

where Fa{r,p) is the phase-space distribution of the parton of type a and I^k) is the collision 
kerne! for the collision terms. The initial state for the parton phase space distribution is taken 
from the expenmantally measured nucleon structure functions and elastic nucleon form-factors, 
as well as nuclear shadowing corrections. The general form of collision term, crkliK){r,p), in the 
above equation includes a subset of all possible n n' interaction proceses. These processes 
consists of 2 -4 2 colisions, 2 -> 1 fusions and 1^2 decays of partons in the lowest order 
perturbative QCD, plus higher order corrections which are associated with multiple emissions 
and absorptions of partons. 

Within this study it observed that isotropy in the momentum distribution is reached very 
quickly, leading to very rapid thermalisation, almost after 0.2/m. of maximum overlap of two 
nuclei at RHIC. The system is very dense {n(r) = 565/m-3(^)-°9°) ^g^y hot {T{r) = 
950 iVfe 1/(5^)“°®°) within this study. 
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8.2 Lattice QCD studies 

III the high temperature limit of QCD, the short range interactions among the partons become 
weak and the long range interactions become dynamically screened. As a consequence, hadronic 
matter at very high temperature do not exhibit confinement and chiral symmetry breaking. 
Lattice gauge theory which allows nonperturbative numerical calculation of observables in QCD 
predicts that the transition between low-temperature and high-temperature phases of QCD is 
not smooth but discontnuous. The order of chiral phase transition is sensitive to the number 
of light, dynamical quark flavors. This phase transition is of second-order for two massless 
flavors and is of first-order for three massless flavors. The numerical simulations of lattice gauge 
theory have predicted that the transition temperature lies in the range 150±10 MeV at vanishing 
net quark density. The simulations without dynamical quarks clearly exhibit a first-order phase 
transition at a temperature Tc ^ 260MeV'. The quarks are not confined in the high temperature 
phase. The quark condensate < > shows a strong drop over the same temperature range, 

indicating that deconfinement and restoration of chiral symmetry go hand-m-hand. However, 
the simulations of lattice QCD with dynamical quarks have not yet overcome the limitations due 
to finite lattice size. More works in this area are now in progress regarding the order of phase 

transition. 
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